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EXECUTIVE  SUMMARY 


The  Installation  Restoration  program  (IRP)  underway  at  numerous  Air  Force 
bases  has  identified  several  sites  with  contaminated  soil  and  groundwater. 

This  subsurface  contamination  is  the  result  of  fuels,  cleaning  solvents,  and 
degreasers  entering  the  subsurface  environment  from  accidental  spills,  leaking 
storage  tanks,  and  past  disposal  practices.  HQ  AFESC/ROVW  is  conducting 
research  aimed  at  developing  treatment  strategies  for  groundwater  cleanup,  and 
studying  the  fate  and  transport  of  contaminants  in  subsurface  systems.  Many 
of  the  contaminants  of  concern  are  volatile  by  nature,  and  a  knowledge  of 
their  air-water  distribution  and  aqueous  solubility  is  needed  to  assess  the 
compounds'  treatability  and  to  support  the  basic  laboratory  studies. 

The  objectives  of  this  research  were  to  develop  Henry's  law  constants  and 
aqueous  solubilities  as  a  function  of  temperature,  for  a  variety  of  organic 
compounds  of  Air  Force  concern  (Table  1).  Secondary  objectives  were  to 
determine  what  effect  mixed  organics,  in  an  aqueous  solution,  exhibit  on 
individual  Henry's  law  constants  and  evaluate  various  methods  used  to  predict 
Henry's  law  constants. 

This  report  documents  experimentally  determined  values  of  Henry's  law 
constants  and  aqueous  solubility  for  51  compounds  of  Air  Force  concern.  The 
report  is  presented  in  three  volumes.  Volume  I  contains  the  technical 
discussion  and  tabulated  values  of  Henry's  law  constants  and  aqueous 
solubilities.  Volumes  II  and  III  contain  all  the  raw  data  and  the  FORTRAN 
source  code  for  an  interactive  program  used  to  predict  the  chemical  parameters. 

Many  of  the  contaminants  of  concern  are  volatile  by  nature,  and  a 
knowledge  of  their  air-water  distribution  is  required  for  the  design  of 
treatment  processes  and  for  providing  insight  into  their  environmental  fate 
and  transport.  A  static  headspace  method  (Equilibrium  Partitioning  In  Closed 
Systems,  referred  to  as  EPICS)  was  used  to  measure  the  Henry's  law  constants, 
with  the  standard  batch  air--stripping  method  used  as  a  check. 

The  Henry's  constants  were  determined  as  a  function  of  temperature  from  10  to 
30  °C  (Tab’e  11)  and  these  values  were  then  used  to  generate  temperature 
regression  equations  (Table  8).  Generally  speaking  the  EPICS'  results  from 
this  study  agree  well  with  other  published  results  (Table  12).  However,  for 
many  of  the  compounds  reported  here,  confirmed  values  of  Henry's  constant  do 
not  exist  in  the  literature,  and,  if  they  do,  values  are  rarely  reported  as  a 
function  of  temperature  with  rigdrous  statistics. 

Solubility  data  for  organic  compounds  in  water  are  important  for 
environmental  studies  because  they  provide  fundamental  information  necessary 
to  predict  transport  in  aqueous  systems.  These  data  may  also  be  used  to 
predict  -arbon  sorption  of  contaminants,  and  the  air-or  steam-stripping 
behavior  for  a  given  compound.  The  aqueous  solubilities  of  the  51  study 


compounds  were  determined  at  10,  20,-  and  30OC  (Table  14).  Three  different 
methods  were  used,  but  the  majority  of  the  data  were  collected  using  a 
shake-flask  technique.  Although  the  solubilities  were  not  a  strong  function 
of  temperature  over  the  range  studied  (I.e.,  10-30  °C),  several  general 
trends  were  noted.  First,  the  solubility  of  the  halongenated  hydrocarbons 
Increased  with  temperature.  Second,  the  solubility  of  the  substituted 
aromatic  hydrocarbons  Increased  with  temperature.  Finally,  maxima  and  minima 
were  observed  for  a  wide  range  of  compounds  without  any  general  trend  that  can 
be  demonstrated  to  be  statistically  significant. 

Groundwater  contamination  Is  often  characterized  by  the  presence  of 
several  different  contaminants,  rather  than  one  single  compound.  For  this 
reason,  studies  were  conducted  to  determine  whether  the  presence  of  other 
compounds  would  affect  the  Henry's  law  constant  of  a  single  compound. 
Deviations  from  Ideal  behavior  were  observed  (pg  52),  but  confirming 
experiments  were  not  performed.  Although  the  results  were  not  conclusive,  the 
project  team  believes  the  observed  Interactions  were  real  and  reproducible. 

It  would  not  be  feasible  to  experimentally  determine  Henry's  law  constants 
for  aP  chemical  compounds.  There  will  be  times  when  a  Henry's  law  constant 
Is  needed  but  an  experimentally  determined  value  Is  not  reported  and  the 
situation  does  not  permit  a  laboratory  study  to  determine  the  constant  For 
this  reason,  a  technique  to  accurately  estimate  Henry's  constant  using  a 
minimum  of  physlochemical  properties  would  be  useful.  Three  different 
thermodynamic  techniques  for  correlating  experimental  Henry's  law  constants 
were  examined  (page  61).  fhe  techniques  were  examined  to  determine  their 
applicability  to  environmental  systems  and  their  predictive  capacity  for 
unmeasured  multicomponent  systems.  The  UNIFAC  method  proved  to  be  the  most 
effective  way  of  utilizing  the  data  base  developed  during  this  project.  A 
co..puter  algorlthim  to  fit  the  current  data  to  a  new  environmental  UNIFAC 
binary  Interaction  data  base  was  developed  and  a  portion  of  the  experimental 
data  collected  was  Incorporated  Into  this  new  data  base.  The  new  data  base 
creates  Improvement  In  the  predictions  generated  by  UNIFAC  In  the  dilute 
concentration  regime  (Figures  13  through  16). 
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The  report  documents  Henry's  law  constants  and  aqueous  solubilities,  as  a 
function  of  temperature,  for  51  compounds  of  Air  Force  concern.  The  study  was 
performed  between  February  1985  and  September  1986. 

This  report  Is  presented  In  three  volumes.  Volume  I  contains  the 
technical  discussion  and  the  tabulated  values  of  Henry's  law  constants  and 
aqueous  solubilities.  Volume  II  contains  the  experimental  Henry's  law  data. 
Volume  III  contains  the  experimental  solubility  data  and  the  FORTRAN  source 
code  for  the  simplex  UNIFAC  parameter  fitting  and  the  Interactive  program  for 
calculating  Henry's  law  constants  and  aqueous  solubilities. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

The  objectives  of  this  research  effort  are  to  develop  Heti'-y's  constants 
and  aqueous  solubilities  as  a  function  of  temperature  for  a  variety  of  organic 
compounds  of  Air  Force  concern.  Secondary  objectives  are  to  determine  what 
effect  mixed  organics  in  an  aqueous  solution  exhibit  on  individual  Henry's  Law 
constants  and  to  correlate  changes  in  chemical  structure  with  Henry's  constants. 
Fifty-one  chemicals  of  Interest  to  the  Air  Force  for  evaluation  of 
air-stripping  technology  were  examined.  The  results  are  presented  in  this 
report  in  considerable  detail,  and  the  methodology  for  development  of  a 
chemical  correlation  is  developed. 

B.  BACKGROUND 

The  current  Installation  Restoration  Program  (IRP)  underway  at  numerous 
Air  Force  bases  has  identified  several  groundwater  contamination  sites.  These 
contaminants  are  the  result  of  fuels,  cleaning  solvents,  and  degreasers  enter¬ 
ing  the  groundwater  from  past  disposal  practices,  accidental  spills,  and  leak¬ 
ing  storage  tanks.  The  Air  Force  is  currently  conducting  an  active  research 
program  to  investigate  feasible  treatment  strategies  for  groundwater  cleanup 
and  to  study  the  interactions  between  groundwater  pollutants  and  soils. 

Packed  tower  air-stripping  is  one  of  the  treatment  strategies  being  devel¬ 
oped  by  the  Air  Force  and  has  proven  to  be  an  excellent  techn’que  for  removing 
trichloroethylene  from  groundwater.  The  same  technique  should  be  applicable 
to  other  volatile  organic  compounds,  but  pilot-scale  testing  is  needed  to 
verify  system  designs.  Accurate  values  of  Henry's  law  constants  are  critical 
in  obtaining  reliable  design  data.  Rarely  will  a  single  compound  be  the  sole 
contaminant  found  in  groundwater.  Groundwater  contamination  is  usually 
characterized  by  the  presence  of  several  of  these  compounds.  For  this  reason, 
studies  are  needed  to  identify  what  effect  complex  mixtures  may  have  on  the 
individual  Henry's  law  constants  for  single  components. 

Henry's  law  constants  have  been  estimated  for  many  compounds  frequently 
using  chemical  property  data  which  are  suspect.  Until  experimentally  deter¬ 
mined  Henry's  constants  are  available,  other  Henry's  law  constants  should  be 
used  with  caution.  In  the  past,  most  Henry's  law  constants  were  experi¬ 
mentally  determined  using  a  bubble-purge  technique  that  is  theoretically 
valid.  However,  the  technique  requires  numerous  assumptions  to  be  practically 
applied.  For  this  reason.  Dr.  J.  Gossett  and  A.-  Lincoff  (Reference  1)  devel¬ 
oped  a  new  method  to  determine  Henry's  law  constants  while  working  on  an  Air 
Force  sponsored  research  project.  The  method,  called  Equilibrium  Partitioning 
In  Closed  Systems  (EPICS),  uses  relative  headspace  concentration  from  two 
equilibrated  systems  to  determine  Henry's  law  constants. 


It  is  not  feasible  to  experimentally  determine  Henry's  law  constants  for 
all  chemical  compounds.  There  will  be  times  when  a  Henry's  law  constant  is 
needed  but  an  experimentally  determined  value  is  not  reported  and  the  situa¬ 
tion  does  not  permit  a  laboratory  study  to  determine  the  constant.  As  men¬ 
tioned  earlier,  techniques  currently  used  to  estimate  Henry's  law  constant, 
such  as  the  ratio  of  vapor  pressure  to  solubility,  have  inherent  limitations  -- 
that  '■estrict  their  use.  For  this  reason,  a  Henry's  law  constant  correlation 
based  on  a  minimum  of  ohysiochemical  properties  would  be  useful.  It  is, 
therefore,  desirable  to  correlate  chemical  structure  with  changes  in  Henry's 
law  constant  for  the  organic  compounds  of  interest  in  this  study,  to  determine 
how  structural  differences  may  influence  equilibrium  partitioning. 

The  Air  Force  is  also  conducting  studies  aimed  at  determining  the  par¬ 
titioning  and  transport  of  volatile  compounds  in  groundwater  systems.  The 
transport  of  contaminants  in  groundwater  systems  is  often  controlled  by  the 
sorption  between  the  contaminant  in  question  and  the  surrounding  aquifier 
materials.  Contaminant  aqueous  solubility  is  an  important  property  used  to 
estimate  sorption  equilibria.  However,  accurate  solubility  data  are  scarce 
and  compounds  are  often  simply  reported  as  soluble,  slightly  soluble,  or 
insoluble.  For  this  reason,  solubilities  determined  as  a  function  of  tem¬ 
perature  for  each  chemical  will  support  sorption  studies.  The  determination 
of  accurate  Henry's  law  constants  will  also  support  these  studies  providing 
reliable  data  which  can  be  used  in  deve’oping  useful  mass  balances  for  chemi¬ 
cal  partitioning. 

C.  SCOPE 

A  comprehensive  literature  survey  was  performed  to  obtrin  Henry's  law 
constants  and  aqueous  solubilities  for  each  chemical  listed  in  Table  1.  The 
literature  values  are  compared  with  the  results  c*  the  research  effort  and 
included  in  this  final  report. 

Henry's  law  constants  were  then  determined  foi  all  the  chemicals  listed  in 
Table  1  using  the  EPICS  technique  when  applicable.  Although  the  EPICS  tech¬ 
nique  has  proven  to  be  a  simple  and  accurate  method  for  Henry's  constant 
determination,  the  EPICS  method  is  reported  to  lose  utility  when  the  dimen¬ 
sionless  Henry's  law  constant  is  greater  than  two  or  three.  Therefore,  the 
bubble  column  method  was  used  to  determine  Henry's  constants  for  several 
chemicals  fcr  comparison  to  the  EPICS  techiique.  The  purpose  of  this  alter¬ 
nate  measurement  was  to  verify  the  measured  Henry's  law  constants. 

Henry's  law  constants  for  each  chemical  listed  in  Table  1  were  determined 
as  a  function  of  temperature  at  10,  15,  20,  25.  and  30  degrees  Celsius.  In 
addition,  the  aqueous  solubilities  for  each  chemical  listed  in  Table  1  were 
determined  as  a  function  of  temperature  at  10,  20,  and  30  degrees  Celsius. 
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TABLE  1.  LIST  OF  STUDY  COMPOUNDS 


Number 

Component  name 

1 

n-Nonane 

2 

n-Hexane 

3 

2-Methylpentane 

4 

Cyclohexane 

5 

1,2-Dichlorobenzene 

5 

Chlo>'obenzene 

7 

1,3-01  chlorobenzene 

8 

1,4-D1 chlorobenzene 

9 

o-Xylene 

10 

p-Xylene 

11 

m-Xylene 

12 

Propyl  benzene 

13 

Ethylbenzene 

14 

Toluene 

15 

Benzene 

16 

Phenol 

17 

Methyl  ethylbenzene 

18 

1,1-Dichloroethane 

19 

1,2-Dichloroethane 

20 

1,1, 1-Trichloroethane 

21 

1,1,2-Trichloroethane 

22 

cis-DIchloroethylene 

23 

trans-Oichloroethylene 

24 

Tetrachloroethylene 

25 

Trichloroethylene 

26 

Naphthalene 

27 

Tetralin 

28 

Decal  in 

29 

Anthracene 

30 

Vinyl  chloride 

31 

Chloroethane 

32 

Hexachloroethane 

33 

Carbon  tetrachloride 

34 

1,3,5-Trlmethylbenzene 

35 

bis(2-Ethylhexyl)  phthalate 

36 

Ethylene  d1 bromide 

37 

1,1-01  chi oroethy 1 ene 

38 

Methylene  chloride 

39 

Chloroform 

40 

1,1,2,2-Tetrachloroethane 

41 

1,2-Dlchloropropane 

42 

Olbromochloromethane 

43 

1, 2, 4-Tr1 chlorobenzene 

44 

2,4-Dlmethylphenol 

45 

1 , 1 , 2-Tr1 chi orotri f 1 uoroethane 

46 

Methyl  ethyl  ketone  (HEX) 

47 

Methyl  isobutyl  ketone  (MIBK) 

48 

Methyl  cellosolve 

49 

bls{2-Chloroethyl)  ether 

50 

Trlchlo-ofluoromethane 

51 

2,6-Dln'trotoluene 
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An  attempt  was  also  made  to  evaluate  the  effects  of  mixed  organics  in 
aqueous  solutions  or,  Individual  Henry's  law  constants.  Based  on  information 
obtained  for  pure  components,  six  compounds,  selected  from  Table  1, 
were  exoected  to  cause  the  greatest  amount  of  deviation  in  individual  Henry  s 
1.'..  onstants.  The  Henry's  law  constants  for  each  compound  in  the  mixture 
wire  determined  and  the  results  were  compared  to  the  pure  component  Henry's 
lav  constants  determined  previously. 

The  final  goal  of  this  study  was  to  derive  a  correlation  to  predict 
Henry's  law  constant  based  on  a  knowledge  of  the  chemical  structure  of  the 
compound  in  question.  Two  approaches  have  been  investigated:  an  extension  of 
the  UNIFAC  solution  of  groups  model  as  proposed  by  Fredenslund,  et  a1.  (Refer¬ 
ence  3);  and  the  parent-derivative  method  of  Leo,  Hansch,  and  Elkins  (Refer¬ 
ence  4).  Once  the  method  is  fully  developed,  the  applicability  of  the  cor¬ 
relation  nay  be  demonstrated  by  selecting  five  chemicals  and  predicting  their 
respective  Henry's  law  constant.  The  predicted  Henry's  law  constants  can  then 
be  compared  to  laboratory  values  determined  earlier.  The  correlation  result¬ 
ing  from  these  efforts  is  relatively  simple  to  apply,  requires  a  minimum  of 
physiochemical  properties,  and  is  adaptable  to  other  hydrophobic  chemicals  and 
mixtures  not  examined  in  tne  current  research  effort. 


SECTION  II 


A  REVIEW  OF  HENRY’S  LAW  CONSTANT  CRITERIA 


Air  and  water  contamination  by  volatile  organic  compounds  (VOCs)  has 
become  a  major  environmental  concern  in  the  United  States,  as  evidenced  by  the 
intensive  monitoring  and  regulatory  activities  of  recent  years.  Current  VOC 
contamination  problems  include  air  emissions  from  hazardous  waste  treatment, 
storage,  and  disposal  facilities  (TSDFs)  and  groundwater  contamination  by 
leaking  underground  fuel  storage  tanks.  In  these  and  other  cases,  a  knowledge 
of  VOC-ai r-water  equilibrium  behavior  is  critical  in  designing  VOC  handling 
and  treatment  processes  that  minimize  environmental  impact.  Because  the  VOCs 
are  usually  present  in  the  liquid  phase  at  very  low  concentrations,  Henry's 
law  for  ideal  dilute  solutions  is  often  appropriate  for  vapor-liquid  equili¬ 
brium  (VLE)  calculations. 

The  ongoing  VLE  program  includes  the  laboratory  determination  of  Henry's 
law  constants  and  aqueous  solubilities  for  hydrophobic  organic  chemicals. 
State-of-the-art  techniques  are  being  employed  for  these  measurements,  and  the 
precision  and  accuracy  of  the  data  are  monitored  by  a  variety  of  statistical 
tests. 

A.  THEORETICAL  SOLUTION  THERMODYNAMICS 

The  governing  equation  for  the  equilibrium  partitioning  of  component  "i" 
between  a  liquid  and  a  vapor  phase  in  contact  is  simply: 

fj  =  f^  (at  equilibrium)  (1) 

where 

f^  =  fugacity  of  component  "i"  in  the  vapor  mixture. 

f^  =  fugacity  of  component  "I"  in  the  liquid  mixture. 

It  IS  not  convenient  to  work  directly  with  fugac’ties,  so  the  following 
expressions  are  usually  introduced  to  relate  physically  measurable  quantities 
in  the  two  phases  at  equilibrium: 

(2) 
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M’.U’.U' 


where 


♦.  ^  fugacity  coefficient  (equal  to  unity  if  the  vapor  behaves  as 
an  ideal  gas) 

7^  =  activity  coefficient  (equal  to  unity  for  an  ideal  solution) 

X|,  =  mole  fractions  of  "1“  in  the  liquid  and  vapor  phases,  respectively 

f°  =  standard  state  fugacity  of  “i". 

Using  Equations  (2)  and  (3)  and  assuming  ideal  vapor  and  liquid  phases,  the 
result  is: 

Y,P=X,f°  .  (4) 

To  apply  this  equation,  the  standard  state  (f?)  must  be  chosen  judiciously, 
based  upon  the  liquid-phase  composition  of  tne  system.  Henry's  law  behavior 
arises  when  the  liquid  mixture  of  interest  is  very  dilute  in  component  "i"  and 
a  direct  proportionality  is  observed  between  the  fugacity  of  “i"  and  its  mole 

fraction  in  the  liquid  phase.  Thus,  f?  becomes  equal  to  the  Henry's  constant. 
Hi,  which  is  defined  by  the  relationship: 


For  Raoult's  Law  behavior  when  the  liquid  mixture  is  almost  pure  component 
“1",  the  standard  state  becomes  the  pure  component  vapor  pressure  (Pj^at).  At 
infinite  dilution,  the  Raoult's  law  relationship  becomes: 

YiP  =  Xi7?PiSat  (6) 

where 

t"  =  the  infinite  dilution  activity  coefficient 
Pi^at  =  the  pure  component  vapor  pressure. 
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Comparing  Equations  (4)  and  (6),  Henry's  law  constant  is  equivalent  to  the 
product  of  the  infinite  dilution  activity  coefficient  and  the  pure  component 
vapor  pressure: 

H,  =  •  (7) 

Hence,  by  knowing  the  pure  component  “i"  vapor  pressure  and  the  Henry's  law 
constant  at  any  given  temperature,  one  may  calculate  the  infinite  dilution 
activity  coefficient  at  that  temperature, 

B.  UNITS  CONVERSION  FOR  HENRY'S  LAW  CONSTANT 

The  Henry's  law  constant  in  Equation  (5)  has  thermodynamic  units  [(atm) 
(kmols  of  liquid)/(kmols  of  gas)]  and  can  be  converted  to  dimensionless  form 
by  multiplying  by  the  ratio  of  the  molar  gas  and  liquid  densities.  This 
second  form  of  Henry's  law  is  useful  for  relating  the  aqueous  concentration  of 
a  hydrophobic  organic  compound  to  its  equilibrium  concentration  in  the  vapor 
phase: 


Cg  =  HjCL  (8) 

where 

Cq  =  molar  concentration  of  organic  in  the  gas  phase,  g  mo1es/m^ 

Cj_  =  molar  concentration  of  organic  in  the  liquid  phase,  g  moles/m^ 

Hj  =  dimensionless  Henry's  law  constant. 

A  third  form  of  Henry's  law  is  often  used,  and  the  Henry's  constant  can  be 
expressed  in  still  another  set  of  units.  The  phase  equilibrium  expression  in 
this  case  relates  the  liquid-phase  concentration  of  "i”  to  Its  partial  pres¬ 
sure  in  the  vapor  phase: 

p,  =  V.p,  =  hJC^  (9) 

where 

=  partial  pressure  of  "1,"  (atm) 

Y.  =  mole  fraction  of  ”1"  in  the  vapor  phase. 

Solving  for  Henry  s  law  constant,  hJ,  g^ves  the  following  equation  and  asso¬ 
ciated  set  of  units: 

HJ  .  !i!l  =  ratm  (m3  of  liquid)]  .  (10) 

‘  ~  I  (kmols  of  gas)  J 

In  summary,  the  required  conversion  equations  between  the  three  different 
sets  of  units  for  Henry's  law  constant  are: 


or 


H,  =  H,(pg)  .  (12) 

All  constants  in  the  text  of  this  report  are  in  the  fortn  of  Equation  (10). 
However,  notice  that  pressure  units  of  both  kilopascals  (kPa)  and  atmospheres 
(atm)  are  employed  for  the  sake  of  scaling  or  comparison  to  other  reported 
data.  The  unit  conversion  between  the  two  is  simply  101.325  kPa/atm.  Appen¬ 
dix  E  list  the  Henry's  law  constants  in  many  sets  of  units  and  in  dimension¬ 
less  form. 

C.  TEMPERATURE  DEPENDENCE  OF  HENRY'S  LAW  CONSTANT 

From  fundamental  thermodynamic  relations,  the  temperature  and  pressure 
dependence  of  Henry's  law  constant  can  be  shown  to  be 


h  -  vr 

d(ln  H,)  =  dT  .  y 


hj  •  molar  enthalpy  of  T  in  the  ideal  gas  state 

*  partial  molar  enthalpy  of  "i"  at  infinite  dilution 

=  partial  molar  volume  of  "1"  at  infinite  dilution. 

For  experiments  conducted  at  constant  pressure,  the  second  term  in  Equation 
(13)  disappears.  The  expression  for  the  temperature-dependency  is  then  ob¬ 
tained  by  performing  an  indefinite  integration  on  the  remainder  of  the  equa¬ 
tion  after  assuming  that  the  "enthalpy  change  of  volatilization"  (h'i  -fi")  is 
constant  with  respect  to  temperature  although  this  assumption  is  not  entirely 
realistic  with  respect  to  some  physical  systems.  The  resulting  equation  is 


Cj  =  integration  constant. 


<,v 


Note  that  the  integration  must  be  done  with  the  H,  or  Hi'  forms  of  Henry's  law 
constant  to  account  for  changes  in  the  bulk  gas  density  with  temperature. 

By  treating  the  quantities  (h'j  -  R“,)/R  and  Cj  as  fitting  parameters, 

Hi  data  obtained  at  various  temperatures  can  be  correlated  to  Equation  (14) 
using  linear  regression  analysis.  The  final  equation  for  H^  has  the  following 
form: 


H 


1 


(15) 


where  A  and  B  are  the  "best-fif  constants  for  the  experimental  data.  All 
compounds  are  currently  being  tested  at  5  ‘C  incements  from  10  *C  to  30  ’C, 
and  a  least-squares  procedure  on  the  resulting  data  is  used  to  determine  A  and 
B.  The  interested  reader  should  refer  to  Appendix  B  of  this  report  for  the 
temperature-regression  plots  generated  during  this  study. 
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SECTION  III 


HENRY’S  LAW  CONSTANT  MEASUREMENTS 


Two  techniques  were  used  to  measure  Henry's  law  constant{s)  for  the  com¬ 
pounds  in  this  study.;  One  was  a  batch  air-strlpping  technique  developed  by 
Mackay  et  al.  (Reference  5).  The  other  was  Equilibrium  Partitioning  in  Closed 
Systems,  a  technique  known  as  EPICS,  originally  proposed  by  Gossett  and 
Lincoff  (Reference  1).  Each  of  these  techniques  will  be  described  in  the 
sections  that  follow,  and  the  results  of  the  measurements  performed  will  be 
summarized. 

A..  BATCH  AIR  STRIPPING 

The  most  common  method  for  laboratory  measurement  of  air-water  partition 
coefficients  for  hydrophobic  organics  is  batch  air  stripping,  first  proposed 
by  Mackay  et  al.  (Reference  5).  Conceptually,  the  procedure  is  simple.  High- 
purity  air  bubbled  through  a  continuously  mixed,  organic-laden  aqueous  solu¬ 
tion  IS  allowed  to  equilibrate  with  the  aqueous  phase.  This  operation  is 
typically  carried  out  in  a  tall  vertical  column  in  which  the  gas  is  Introduced 
at  the  bottom  through  a  fine-pore  frit.  The  combination  of  a  tall  column  and 
a  small  bubble  size  (low  bubble  buoyancy)  increase  the  degree  to  which  gas- 
liquid  equilibrium  is  approached.  By  monitoring  the  equilibrium  gas  (or 
liquid)  concentration  versus  time,  Henry's  law  constant  can  be  determined  with 
a  high  degree  of  precision. 

1.  Governing  Theory 

Henry's  law  constant  is  calculated  by  correlating  concentration  data 
against  a  solute  mass-balance  equation  for  the  column.  The  differential  form 
of  the  mass  balance  is  written: 


IN  -  OUT  =  ACCUMULATION 


where 

Vo 

Y(t) 

G 

C{t) 

V 


=  organic  concentration  at  gas  inlet,  g  moles/m^ 

=  time-variable  outlet  gas  concentration,  g  moles/ra^ 
=  volumetric  gas  flow  rate,  m^/minute 
=  time-variable  liquid  concentration,  g  moles/m^ 

=  total  volume  of  batch  liquid  in  column,  m^. 


(16) 
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By  substituting  the  Henry's  law  relationship  Y{t)  =  h,C(t)  into  Equation  (16) 
and  integrating,  the  following  expression  results: 

H  G 

InC(t)  =  InC^  +  -i-  t  (17) 

where 

t  -  elapsed  time  from  start  of  the  run  (minutes). 

According  to  Equation  (17),  a  logarithmic  plot  of  liquid-phase  (or  gas 
phase)  concentration  data  as  a  function  of  elapsed  time  will  yield  a  straight 
line  with  a  slope  directly  related  to  Henry's  law  constant.  Linear  regression 
analysis  of  the  raw  data  Is  typically  employed  to  determine  the  “best-fit" 
slope.  It  IS  not  necessary  to  know  the  Initial  batch  concentration,  Cq,  be¬ 
cause  only  the  y-intercept  of  the  plot  Is  influenced.  The  preparation  of 
quantitative  calibration  standards  to  precisely  determine  absolute  concen¬ 
trations  is  not  necessary  for  the  same  reason.  Therefore,  any  proportional 
measure  of  liquid  concentration,  such  as  Ultraviolet  (UV)  absorbance  in  the 
water  or  gas  chromatograph  peak  area  from  the  air,  can  be  Inserted  for  C(t)  In 
Equation  (17).  For  the  special  case  in  which  llould  sampling  significantly 
affects  the  column  volume.  Equation  (15)  may  be  written  for  discrete  sampling 
Intervals: 

EAt, 

7^  (18) 

where 

Cj  =  liquid  concentration  at  the  end  of  the  i^  time  interval, 
gmoles/m8 

Atj  =  elapsed  time  on  the  1^  sampling  period.  In  minutes 

V.  =  column  liquid  volume  during  the  1^  sampling  period,  m^. 

This  equation  Is  appropriate  for  measurements  of  liquid  samples  periodically 
withdrawn  from  the  column  during  operation. 

2.  Laboratory  Procedures 

A  bubble-purge  column  was  assembled  from  two  glass,  water- jacketed 
columns,  each  2  feet  long  with  a  l-1nch  inside  diameter.  The  columns  were 
Joined  with  a  TefloiP  union  and  rubber  o-ring  seals.-  A  glass-fritted  disc 
with  a  25  to  50-micron  pore  size  was  press-fitted  into  the  column  Inlet.  Air 
was  supplied  to  the  column  from  a  tank  of  high  purity  compressed  air  with  the 
flow  rate  controlled  by  a  rotameter  and  integral  needle  valve.  Purge  air  was 
sparged  through  a  water-filled  impinger  before  entering  the  bubble  column  to 
prevent  water  evaporation  from  the  column.  The  column  exit  was  fitted  with  a 
glass  tee  for  exhausting  column  effluent  to  the  fume  hood  and  for  sampling  the 
exit  gas  by  syringe  from  a  gas-tight  Teflon-faced  septum.  A  constant  flow  of 
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temperature-controlled  water  was  maintained  through  the  column  jacket  to 
ensure  isothermal  operation  over  the  entire  liquid  depth.  Figure  1  shows  a 
schematic  of  this  system. 


Before  operation,  the  column  was  loaded  with  a  known  volume  of  dis¬ 
til  led/deionized  water.  After  the  water  reached  the  temperature  of  the  water 
jacket,  the  liquid  volume  was  incieased  to  Its  ultimate  value  by  pipetting  a 
quantity  of  saturated  organic  solution  into  the  column.  The  organic  solution 
was  introduced  subsurface  to  minimize  volatilization  losses.  The  use  of  satu¬ 
rated  stock  solution,  rather  than  neat  organic,  eliminates  the  difficulty  of 
getting  the  organic  Into  solution  before  starting  a  run.  Gas  flow  was  then 
initiated  at  the  desired  rate,  and  a  few  minutes  were  allowed  for  column  flow 
patterns  to  stabilize  before  sampling  was  begun. 

Samples  of  column  effluent  gas  were  collected  with  a  gas-tight  syringe 
and  injected  into  a  gas  chromatograph  with  a  flame  ionization  detector.  The 
GC  column  consisted  of  a  14-inch  by  l/8-inch  O.D.  stainless  steel  tube  packed 
with  80/100  mesh  PorasiiP  B.  The  advantage  of  using  this  column  was  that 
sample  peak  elution  occurred  within  30-45  seconds  of  Injection,  allowing  real¬ 
time  analysis  of  column  exit  gas  concentrations. 

For  those  compounds  exhibiting  fast  decay  of  gas-phase  concentration, 
the  variation  in  time  between  sample  collection  and  injection  by  syringe 
significantly  degraded  the  test  reproducibility  and  decay-curve  linearity. 

For  these  compounds,  the  column  exit  gas  was  pumped  through  a  1/16-inch  outer 
diameter  stainless  steel  tube  and  passed  through  a  0.5  cc  stainless  steel 
sample  loop  attached  to  a  six-port  rotary  valve.  This  valve  was  pneumatically 
activated  to  load  the  loop  contents  onto  the  GC  column.  Such  an  arrangement 
allowed  almost  instantaneous  injection  of  sampled  effluent  gas. 

3.  Results  and  Discussion 

Initial  tests  performed  with  the  bubble  purge  column  were  aimed  at 
optimizing  the  operating  conditions  and  experimental  procedures.  Hexane  was 
chosen  as  the  first  test  compound  because  of  its  high  Henry's  law  constant,  as 
calculated  from  vapor  pressure  and  solubility  data.  Hexane's  high  volatility 
would  likely  reveal  any  problems  in  reaching  equilibrium  with  the  chosen 
column  design. 

Tests  were  performed  at  a  column  temperature  of  30  ‘C  with  the  flow 
rate  being  varied  between  tests.  The  results  of  several  tests  performed  using 
hexane  solutions  are  shown  in  Table  2. 


TABLE  2.  BUBBLE  PURGE  TEST  RESULTS— HEXANE 


Air  Flow  Rate,  cc/mm 

Liquid  Depth,  cm 

He,  kPa-m3/mol 

91.6 

124 

5.948 

47.6 

121 

8.014 

48.4 

119 

7.608 

21.8 

123 

14.913 

14.1 

128 

19.050 

12 


AS 

At 

•Cr 


At  the  lower  gas  flow  rates,  a  smaller  bubble  size  was  produced  and  there 
appeared  to  be  a  corresponding  improvement  in  the  approach  to  equilibrium  as 
evidenced  by  the  increasing  Henry's  law  constant.-  However,  a  subsequent  test 
with  only  one  section  of  the  bubble  purge  column  in  place  and  a  liquid  depth 
OT  5b  cm  produced  a  calculated  Henry's  law  constant  of  7.540  kPa-m^/mole, 
nearly  as  high  as  the  value  observed  at  a  depth  of  119  cm  with  the  same  purge 
gas  flow  rate.  If  the  changes  observed  in  the  Henry's  law  constant  with 
changing  gas  flow  rate  were  caused  by  equilibrium  conditions  not  being  met, 
one  would  expect  a  significant  effect  on  the  constant  from  decreasing  the 
liquid  depth  by  one-half. 

Further  tests  were  performed  with  1, 1 , l-trich1oroethane  solutions. 
This  compound  has  a  much  lower  calculated  constant  than  hexane  and  should  have 
been  much  less  susceptible  to  column  equilibrium  limitations.  Tests  were  all 
performed  at  a  column  temperature  of  30  *C  with  the  liquid  depth  and  gas  flow 
rate  being  varied  between  tests.  These  results  are  shown  in  Table  3. 


TABLE  3.  BUBBLE-PURGE  TEST  RESULTS-1,1,1-TRICHL0R0ETHANE 


Air  Flow  Rate,  cc/rain 

Liquid  Depth,  cm 

H,  kPa-m^/mol 

24.8 

54 

2.258 

27.2 

120 

3.636 

27.3 

122 

3.508 

104.0 

54 

2.212 

106.0 

75 

2.474 

106.0 

100 

2.756 

108.0 

120 

2.870 

Two  observations  can  be  made  from  these  results.  First,  for  a  par¬ 
ticular  flow  rate.  Increasing  the  liquid  depth  Increased  the  Henry's  law  con¬ 
stant  measured.  Second,  a  variation  in  the  measured  Henry's  law  constant  with 
changes  in  flow  rate  at  a  particular  liquid  depth  was  observed  only  for  tne 
two-section  column  (54  cm  liquid  depth  tests  were  conducted  with  a  single- 
column  section).  This  second  axial  mixing  within  the  column  resulted  in  a 
bias  in  the  measured  decay  rate,  because  such  an  effect  would  be  more  pro¬ 
nounced  at  lower  flow  rates  and  greater  liquid  depths. 

To  find  column  mixing  efficiency,  a  test  was  performed  to  observe  the 
dispersion  of  food  coloring  with  the  column  liquid  during  operation.  Water 
(aoout  575  mL)  was  added  to  the  column  to  a  depth  of  115  cm  and  air  flow  was 
established  at  nominally  25  cc  per  minute.  A  few  drops  of  red  food  coloring 
were  added  to  the  water  at  the  top  of  the  column  and  the  dispersion  was 
observed.  Complete  mixing  was  observed  in  the  top  column  section  within  2 
minutes;  however,  approximately  25  minutes  elapsed  before  the  bottom  column 
section  reached  the  same  color  intensity. 

Two  steps  were  taken  to  prevent  axial  mixing  problems  from  biasing 
the  results  for  the  remaining  measurements  using  batch  air  stripping.  First, 
only  the  lower  2-foot  column  section  was  used  because  of  the  considerable 
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TABLE  6.  BATCH  AIR-STRIPPING  RESULTS— HEXANE 


Temperature,  *C 

Liquid  depth,  cm 

H,  kPa-ml/mo1 

10 

54.0 

4  468  +  0.054  (1.2) 

10 

57.5 

4.074  +  0.016  (0.4) 

15 

54.0 

4.700  +  0.141  (3.0) 

15 

57.5 

4.780  +  0.303  (6.3) 

20 

54.0 

5.265  +  NA 

20 

57.5 

5.168  +  0.026  (0.5) 

25 

54.0 

5.391  +  0.032  (0.6) 

25 

57.5 

5.397  +  0.167  (3.1) 

30 

54.0 

5  957  f  0.038  (0.6) 

30 

57.5 

5.973  +  NA 

TABLE  7. 

BATCH  AIR  STRIPPING 

RESULTS-METHYL  PENTANE 

Temperature,  *C 

Liquid  depth,  cm 

H,  kPa-m^/mol 

10 

54.0 

4.322  +  0.134  (3.1) 

10 

57.5 

4.282  +  0.076  (1.8) 

15 

54.0 

4.618  +  0.098  (2.1) 

15 

57.5 

4  740  +  0.040  (0.85) 

20 

54.0 

5.018  +  0.052  (l.O) 

20 

57.5 

5.010  +  0.099  (2.0) 

25 

54.0 

5.322  0.014  (0.27) 

25 

57.5 

5.292  +  0.176  (3.3) 

30 

54.0 

5.990  +  0.004  (0.06) 

30 

57.5 

5.956  +  0.009  (0.15) 

B.  EQUILIBRIUM  PARTITIONING  IN  CLOSED  SYSTEMS  (EPICS) 


1,'  Governing  Theory 


Because  of  the  equilibrium  limitations  of  the  batch  air-stripping 
technique,  the  EPICS  method  is  an  attractive  alternative  for  determining 
Henry's  law  constant.  The  EPICS  procedure  (Referenc;  1)  is  based  on  a  closed 
system  mass  balance  for  a  given  VOC  distributed  between  liquid  and  gas  phases 
in  contact.  A  component  balance  takes  the  general  form: 


where 


M  =  C  CgVg 


(19) 


M  =  total  organic  mass  in  the  system 
Cj^  =  liquid-phase  organic  concentration 


Cg  =  gas-phase  organic  concentration 

V|^  =  total  liquid  volume 

Vg  =  total  gas  volume 

If  the  same  mass  of  organic  is  introduced  into  two  closed  containers 
(i.e.,  sealed  septum  bottles)  containing  different  volumes  of  pure  water,  an 
expression  similar  to  Equation  (19)  may  be  written  for  eacn  system.  Equating 
these  mass  balance  expressions  and  introducing  Henry's  law  to  substitute  for 
the  liquid-phase  VOC  concentrations  gives  an  equation  relating  Henry's 
constant  to  headspace  concentrations  and  known  volumes: 


H' 


/Cg  )V  -  V 
'’2  '-1  “-2 

-  (C.  /Cr  )V. 

Gj  Gj 


(20) 


where 


H'  =  dimensionless  Henry's  law  constant,  atm  •  (m^  of  liquid)/(m3 
of  gas) 

1,2  =  subscripts  identifying  the  two  closed  systems. 

The  only  experimental  information  needed  to  determine  Henry's  law  constant  is 
the  ratio  of  gas-phase  concentrations  (C-  /C.  )  in  the  two  systems.  In  fact, 

bj  b^ 

absolute  concentrations  are  not  necessary  because  any  proportional  measure  of 
concentration  (such  as  GC  peak  areas)  will  yield  the  desired  headspace  ratio 
in  Equation  (20)., 

2.  Laboratory  Procedures 

In  practice,  the  EPICS  method  involves  measurement  of  multiple  pairs 
of  high/low  volume  closed  systems  to  obtain  several  independent  estimates  of 
Henry's  law  constant.-  If  three  pairs  are  used,  for  example,  three  replicate 
headspace  measurements  will  be  obtained  for  each  of  the  two  liquid  volumes. 
Thus,  the  six  values  can  be  paired  in  all  permutations  to  calculate  nine 
estimates  of  Henry's  law  constant.  Typically,  an  arithmetic  average  of  these 
estimates  is  taken  to  be  the  true  Henry's  law  constant,  with  the  coefficient 
of  variation  for  the  individual  values  reported. 

The  EPICS  technique  used  during  the  current  VLE  program  uses  two 
replicate  bottle  pairs  per  observation  to  obtain  four  estimates  of  Henry's 
constant,  and  a  second  set  ot  selected  compounds  is  analyzed  on  a  separate  day 
to  deteiTBine  day-to-day  precision.  Therefore,  a  total  of  eight  Henry's  law 
constant  estimates  are  obtained  for  a  selected  component. 

A  saturated  stock  solution  fcr  each  component  is  prepared  by  adding  an 
amount  of  organic  solute  slightly  in  excess  of  the  solubility  limit  to  a 
volume  of  pure  deionized/disti Ilea  water.  AH  organic  chemicals  were  used  at 
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least  99  percent  purity.  All  stock  solutions  are  prepared  in  1-liter  amber 
bottles  and  allowed  to  equilibrate  for  a  minimum  of  1  week  before  use.  For 
compounds  more  dense  than  water,  solution  is  withdrawn  directly  from  the  amber 
bottles  to  prepare  EPICS  samples.  For  compounds  less  dense  than  water,  the 
stock  solution  is  transferred  to  a  500  ml  separatory  funnel  the  day  before  the 
scheduled  EPICS  tests,  then  withdrawn  the  next  day  from  the  bottom  of  the 
separatory  funnel  to  prepare  the  EPICS  sample  bottles. 

In  preparing  EPICS  samples,  four  250  mL  amber  glass  bottles  are  filled 
with  20  ml  and  200  ml,  respectively,  of  deionized/distllled  water  into  two 
bottles.  The  same  volume  of  saturated  stock  solution  is  then  added  to  each  of 
the  four  septum  bottles.  The  exact  solution  volume  added  is  dependent  upon 
compound  so'ubility  and  each  test  was  chosen  to  produce  an  initial  liquid 
concentration  of  10  mg  per  liter  in  the  200  ml  system.  For  high  concentration 
stock  solutions,  1  cc  was  chosen  as  the  lower  limit  for  addition  to  the  EPICS 
bottles..  For  low  concentration  stock  solutions,  10  cc  was  chosen  as  the  upper 
limit  to  be  added. 

In  addition  to  the  EPICS  san^les,  one  extra  bottle  per  compound  Is 
prepared  for  use  in  testing  and  conditioning  the  GC  system.  Also,  one  blank 
bottle  consisting  of  deionized/distllled  water  only,  Is  prepared  for  everv  10 
EPICS  sample  bottles  prepared.  The  bottles  are  sealed  with  a  silicone  rubber 
spectrum  cap  with  a  Teflon®  liner  facing  toward  the  bottle  headspace.  These 
liners  are  used  only  once  to  prevent  adsorption  of  test  compound  into  the 
silicone  rubber  once  the  liner  has  been  pierced.  The  loaded  bottles  are 
shaken  vigorously  by  hand  and  then  completely  submersed  in  a  constant  tempera¬ 
ture  water  bath  for  a  minimum  of  16  hours  before  analysis.  This  period  was 
determined  to  be  adequate  to  attain  a  constant  (equilibrium)  headspace  concen¬ 
tration  with  respect  to  time  for  all  of  the  chemicals  examined. 

After  equilibration,  headspace  samples  are  withdrawn  from  the  bottles 
via  Hamilton  1.0  cc  gas-tight  syringe  and  injected  into  a  Varlan  3700  gas 
chromatograph  equipped  with  a  flame  ionization  detector  (FID).  The  GC  column 
is  Identical  to  that  used  for  batch  air-stripping  tests.  Component  retention 
times  are  typically  less  than  1  minute. 

To  check  detector  response  linearity  and  the  validity  of  Henry's  law 
over  the  concentration  range  used,  a  series  of  quantitative  standards  is 
prepared  and  analyzed  for  each  component.  The  standard  series  consists  of 
varying  volumes  of  saturated  stock  solution  diluted  to  the  same  volume  with 
deionized/distllled  water  in  250  ml  septum  bottles.  Response  linearity  is 
checked  by  plotting  headspace  gas  chromatograph  (GC)  response  versus  volume  of 
saturated  stock  solution.  The  results  of  all  of  these  tests  are  shown  in 
Appendix  B. 

3.  Results 

Measurement  of  Henry's  law  constant  using  EPICS  has  been  performed  for 
48  of  the  51  compounds  proposed.  The  headspace  concentrations  produced  by 
samples  of  anthracene,  bis  (2-ethylhexyl )  phthalate,  and  1,2-dinitrotoluene 
were  below  the  limit  of  detection  of  the  GC-FIO.  This  was  expected,  based  on 


the  extremely  low  solubility  (ana  thus,  low  stock  solution  concentration)  of 
these  compounds  as  well  as  low  vapor  pressure.  Unfortunately,  liquid-phase 
measurements  for  these  materials  were  also  inadequate,  because  virtually  no 
organic  partitions  into  the  headspace,  and  the  measurements  cannot  be  distin¬ 
guished  from  the  dilution  factors.  For  the  remaining  compounds,  four  con¬ 
stants  were  determined  for  each  compound  at  temperatures  of  10,  15,  20,  25, 
and  30  *C.  These  results  are  summarized  in  Table  8  in  terms  of  the  tempera¬ 
ture  regression  study  and  are  presented  in  detail  in  Appendix  8. 

Measurements  of  Henry's  law  constant  have  been  repeated  for  these 
compounds  that  exhibited  poor  correlation  coefficients  for  the  linear  regres¬ 
sion  of  H,  (Henry's  law  constant)  versus  ■•eciprocal  absolute  temperature.  In 
the  repeat  measurements,  EPICS  bottles  were  immersed  in  the  water  bath  to 
their  necks  during  equilibration.  This  was  made  possible  by  preparing  fewer 
bottles  and  using  metal  clips  to  prevent  low  liquid  volume  bottles  from 
floating. 

A  major  objective  of  this  investigation  was  to  reduce  the  Henry's  law 
constant  data  to  a  final  form  useful  for  manual  or  computer  calculations.  The 
chosen  means  to  accomplish  this  was  to  generate  temperature  regression  equa¬ 
tions  based  on  the  theoretical  constant-pressure  temperature  dependence  of 
Henry's  law  constant.  The  expression  describing  this  dependence  was  presented 
earlier  as  Equation  (15).  The  “enthalpy  of  volatilization"  is  assumed  con¬ 
stant  over  the  entire  temperature  range  in  the  regression,  even  though  such  a 
simplification  may  not  always  be  physically  realistic.  Generally,  the  simple 
two-constant  fit  seemed  to  represent  the  data  well,  as  evidenced  by  typical 
r-squared  values  in  excess  of  0.95.  Given  the  rigorous  statistics  calculated 
for  the  data  over  the  course  of  the  study,  attaching  reliable  error  bounds  to 
the  regression  equation  proved  possible.  Table  8,  therefore,  completely 
describes  the  data  as  a  temperature  regression  equation  coupled  with  a  tem¬ 
perature-dependent  error  term,  with  upper  and  lower  error  limits  in  absolute 
units  (Henry's  law  constant  units)  or  as  upper  and  lower  percentage  limits. 

Of  the  original  51  chemicals,  the  results  for  48  are  presented  in  Table  8. 

An  examination  of  Table  8  reveals  very  good  r-squared  values  for  the 
linear  temperature  regression  for  most  of  the  chemicals  examined;  however, 
such  results  have  not  been  achieved  for  a  few  compounds.  Using  an  r-squared 
value  of  0.95  as  a  somewhat  arbitrary  cutoff  point  for  assessing  the  quality 
of  it,  the  initial  temperature  regression  fits  for  21  of  the  51  chemicals  of 
interest  were  deemed  to  be  inadequate.  The  entire  EPICS  series  was  repeated 
for  these  21  chemicals.  Careful  replicatio.i  of  these  compounds  in  some  cases 
yielded  correlations  of  improved  quality.  This  can  be  taken  as  an  indication 
that  the  EPICS  technique  itself  is  not  the  cause  of  the  poor  temperature 
regression  correlation  noted  for  a  few  of  the  chemicals  Investigated,  but  that 
either  a  linear  regression  equation  does  not  appropriately  describe  the 
thermodynamics  for  that  mixture,  or  the  Henry's  law  constant  is  outside  the 
range  appropriate  to  headspace  measurements.  Iinprovement  for  some  of  the 
compounds  can  be  attributed  to  the  refinement  of  the  E^ICS  technique  over  the 
project's  course.  Table  9  shows  the  95  percent  confidence  bands  for  the  48 
chemicals  25*  C  calculated  from  the  equations  presented  in  Table  8. 
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TABLE  8.  TEMPERATURE  REGRESSION  PARAMETERS  AMD  ERROR  BAM>S  (AT  26  *C)  FOR  61  ORGANIC  CHEMICALS 
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In  some  cases,  the  linear  correlation  of  the  logarithm  of  the  Henry's 
law  constant  and  temperature  as  described  by  Equation  (13)  results  In  an 
exceedingly  poor  correlation  coefficient,  even  for  replicated  compounds.  For 
example.  In  Figure  2  the  tempcature  regression  plot  for  2-methyl  pentane 
(Component  103)  is  shown.  For  this  compound  a  linear  regression  may  not  be 
appropriate,  but  clearly  the  data  could  be  represented  by  a  smooth  curve  with 
a  local  minimum  near  20  *€. 

The  Henry'S  law  constant  may  be  defined  In  terms  of  the  infinite 
dilution  activity  coefficient  and  t.he  pure-component  vapor  pressure  (Equation 
7).  Many  solutions,  particularly  aromatics  in  water,  are  known  to  have  an 
activity  coefficient  maximum  in  the  range  of  10-25  *C  C^sonopaulos  and 
Praunsnitz,  Reference  6).  This  being  the  case.  It  is  not  surprising  that 
those  materials  with  a  strong  activity  coefficient  dependence  on  temperature, 
might  also  show  a  strongly  nonlinear  Henry's  law  constant  behavior.  As  shown 
in  Section  III  on  the  solubility  measurements,  minima  and  maxima  in  solubility 
are  indeed  noted  with  increasing  temperature.  A  theoretically  sound  correla¬ 
tion  of  activity  coefficients,  taking  into  account  the  structure  of  water  as  a 
function  of  temperature,  has  not  yet  been  proposed  to  explain  this  phenomenon. 

A  second  and  probably  more  important  measure  of  data  quality  is  the 
coefficient  of  variation  (COV)  for  individual  data  sets,  particularly  for 
replicates  at  one  temperature.  However,  poor  temperature  regressions 
generally  correspond  to  high  COV  values  for  the  EPICS  four-bottle  series  at 
all  temperatures.  Figures  3  through  7  present  the  COV  values  for  the 
compounds  with  a  good  linear  regression  correlation  plotted  against  their 
respective  Henry's  law  constants  at  each  of  the  five  target  temperatures  (10, 
15,  20,  25,  and  30  *C).  Table  10  summarizes  all  of  the  COV  values.  The  only 
discernable  trend  apparent  in  these  plots  is  a  general  rise  in  the  COV  value 
with  decreasing  volatility,  an  effect  that  may  be  due  to  reduced  sensitivity 
of  the  analytical  Instrumentation.  An  examination  of  the  average  COV  values 
In  Figures  3  through  7  Indicates  an  overall  EPICS  four-bottle  precision  of  5 
percent  or  better  over  the  whole  volatility  range. 

Finally,  the  actual  Henry's  law  constants,  as  measured  at  each  of  the 
temperatures  for  all  of  the  chemicals,  are  listed  in  Table  11.  The  expanded 
raw  data  including  all  of  the  statistical  data  are  listed  In  Appendix  B. 

4.  Evaluation  of  the  Results 

A  comparison  of  selected  results  of  this  study  and  results  of  previous 
studies  is  found  in  Table  12.  Comparisons  of  Henry's  law  constants  from  EPICS 
'n  this  study  to  bubble  column  results  are  quite  good  with  three  exceptions. 

A  difference  of  over  69  percent  is  observed  for  cyclohexane  and  can  be  attri¬ 
buted  to  the  failure  to  attain  equilibrium  In  the  bubble  column  tests. 

As  expected,  the  greatest  difference  between  the  results  of  this  study 
and  other  reported  values  occurs  In  the  comparison  with  Henry's  Law  constants 
predicted  from  vapor  pressure  and  solubility  data.  The  differences,  both 
positive  and  negative,  ranged  from  2  to  400  percent,  with  no  distinct  pattern 
attributable  to  chemical  type  or  structure.  However,  It  is  not  uncommon  for 
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Figure  2.  Temperature  Regression  Plot  for  2-Methyl  Pentane. 


EPICS  lochnique.  10  C 


Figure  3.  Coefficient  of  Variance  Values  v$.  Henry's  Constant  for  the  EPICS  at  10 
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Figure  5.  Coefficient  of  Variance  Values  v$.  Henry's  Constant  for  the  EPICS  at  20 
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reported  values  of  Henry's  constants  to  vary  by  several  orders  of  magnitude. 
For  example,  experimental  results  for  vinyl  chloride  suggest  a  Henry's  law 
constant  of  approximately  2.21  by  10‘2  atm  by  m^/mol  at  25  *C,  while  values  of 
2.3  by  10"2  to  6.39  atm  by  m^/mol  are  reported  by  Mackay  and  Shiu 
(Reference  2)  and  Goldstein  (Reference  11).  Large  uncertainties  in  the 
absolute  aqueous  solubility  probably  contribute  most  to  these  discrepancies.. 

The  comparison  of  the  EPICS  results  from  this  study  to  those  of 
Gossett  (Reference  7)  are  '"ilte  good.  Again,  both  positive  and  negative 
differences  were  observ  all  values  agreed  within  10  percent.  This 

close  agreement  demonstra.'  cw  reproducible  results  are  between  investi¬ 
gators  using  the  EPICS  tech.i  ^ue. 

Leighton  and  Calo  (Reference  8)  used  an  equilibrium  cell  with  con¬ 
tinuous  gas  flow,  followed  by  direct  analysis  of  both  the  air  and  water  phases 
to  determine  Henry's  constants.  Of  the  nine  compounds  common  to  both  studies, 
results  for  six  of  the  compounds  agree  reasonably  well  (within  10  percent). 
Larger  differences  are  observed  for  the  other  three  compounds  (1,1,1- 
trichloroethane,  1,1,2-trlchloroethane,  and  1,1-dichloroethylene).  Both 
positive  and  negative  differences  (15  to  42  percent)  were  noticed.  This 
difference  can  be  explained  by  the  analytical  problems  discussed  in  Leighton's 
study  and  is  not  an  artifact  of  the  EPICS  procedure. 

The  UNIFAC  (UNIFAC  Functional  Group  Activity  Coefficient)  model  was 
developed  in  1975  to  correljte  large  quantities  of  data  with  a  relatively  few 
molecular  parameters.  The  first  group  contribution  technique  developed  for 
activity  coefficient  prediction  was  that  of  Wilson  and  Deal  (Reference  12); 
the  subsequent  development  of  UNIFAC  by  Fredenslund  et  al.  (Reference  3)  owes 
much  to  their  early  work.  The  theoretical  aspects  of  UNIFAC  are  discussed  in 
detail  in  Section  VI. 

By  employing  the  UNIFAC  model  to  calculate  activity  coefficients  for 
the  compounds  of  interest  in  a  dilute  aqueous  solution,  one  may  then  use 
Equation  (7)  in  conjunction  with  pure  component  vapor  pressure  data  to  calcu¬ 
late  Henry's  law  constants.  Values  obtained  in  this  fashion  are  shown  in  the 
last  column  of  Table  12.  A  limitation  of  this  approach  is  that  the  binary 
intei’action  paraincters  tabulated  by  Gmehling  et  al,  (Reference  13)  were  not 
derived  from  data  in  the  extremely  dilute  region.  As  a  result,  extrapolation 
to  infinite  dilute  may  lead  to  large  errors,  as  seen  for  cyclohexane  or  tetra- 
chloroethylene. 

Generally  speaking,  the  EPICS  results  f'cim  this  study  agree  well  with 
other  published  results.  However,  for  most  of  the  compounds  reported  here, 
reliable  values  of  Henry's  constant  do  not  exist  in  the  literature,  and  if 
they  do,  ‘/alues  are  rarely  reported  as  a  function  of  temperature  with  rigorous 
statistics.  The  EPICS  procedure  proved  to  be  a  sinple  and  accurate  technique 
for  determining  Henry's  law  constants  for  a  variety  of  organic  chemicals.  The 
constants  reported  have  numerous  environmental  applications  (fate  and  trans¬ 
port  models,  remedial  action  plans)  and  are  within  the  10  percent  accuracy 
suggested  by  Mackay  and  Shiu  (Reference  2). 
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SECTION  IV 


SOLUBILITY  MEASUREMENTS 


A.  INTRODUCTION 

Solubility  data  for  organic  compounas  in  water  are  important  for  envi¬ 
ronmental  studies  because  they  provide  fundamental  information  necessary  to 
predict  transport  in  aqueous  systems.  This  data  may  also  be  used  to  predict 
carbon  adsorption  of  contaminants,  and  the  air  or  steam  stripping  behavior  for 
a  given  organic.  Tor  highly  soluble  materials  these  data  are  readily  obtained 
by  methods  as  simple  as  weight  measurements;  however,  precise  determination  of 
solubility  limits  is  difficult  for  sparingly  soluble  compounds,  requiring 
tedious  methods  such  as  radio  labelling  {Burns  and  MacIntyre,  Reference  14} 
to  obtain  precise  results.  Rarely  is  any  temperature  dependency  data  avail¬ 
able  for  solubility;  therefore,  in  this  report,  aqueous  solubility  of  the  5! 
tests  chemicals  is  reported  at  10  *C,  20  *C,  and  30  *C. 

Three  methods  dominate  the  current  literature--nephelometry,  the  shake- 
flask  method,  and  the  generator  column— and  each  has  specific  limitations.- 
Nephelometry  differs  from  the  latter  two  in  that  the  measurement  employed  is 
indirect.  Mixtures  of  water  and  the  solute  of  interest  are  made  for  several 
known  dilutions.  These  mixtures  are  agitated  vigorously,  using  either 
mechanical  shakers  or  ultrasonic  devices  to  produce  a  total  suspension  of  the 
solute  in  water.  The  turbidity  of  these  mixtures  is  then  determined  optically 
for  a  succession  of  solute  dilutions  and  plotted  as  a  function  of  the  known 
concentration  of  these  dilutions.  The  line  obtained  is  extrapolated  to  zero 
turbidity.  The  concentration  at  that  point  is  taken  as  the  solubility  limit 
of  the  solute.  Hollifield  (Reference  15)  reports  95  percent  confidence  limits 
of  between  5  and  20  percent,  with  decreasing  accuracy  for  less  soluble 
materials,  using  this  method.  However,  the  employed  surfactants  and  other 
materials  to  achieve  uniform  emulsions.  Because  of  the  poor  accuracy  for 
highly  insoluble  chemicals,  most  researchers  have  chosen  to  use  the  shake- 
flask  or  column  generator  method  to  determine  the  solubilities  of  these 
compounds. 

In  the  shake-flask  method,  a  known  amount  of  solute  is  added  to  a  measured 
water  sample  in  a  sealed  container.  The  amount  of  solute  is  chosen  to  exceed 
the  anticipated  solubility  limit.  This  mixture  is  shaken  or  stirred  and 
allowed  to  sit  for  3  hours  to  several  days  to  ensure  that  solute-solvent 
equilibrium  is  obtained.  The  mixture  is  then  filtered  or  centrifuged  to 
remove  suspended  solute  particles,  and  the  concentration  in  the  remaining 
solution  is  determined. 

Two  basic  methods  may  be  used  to  determine  the  resulting  aqueous  concen¬ 
tration.  One  method  is  to  employ  Henry's  law  and  the  measurement  of  the  gas- 
phase  concentrations  in  the  headspace  to  determine  the  aqueous  concentration. 


The  gas-phase  concentration  of  the  solute  is  determined  for  the  saturated 
solution  and  divided  by  the  Henry's  law  constant  to  obtain  the  aqueous 
solubility  limit. 

The  second  method  of  determining  the  aqueous  concentration  involves 
extracting  the  dissolved  material  from  a  known  volume  of  water  and  determining 
the  amount  of  materials  recovered.  The  water  solution  is  passed  through  a 
sorbent  material  bed  to  remove  any  organ'c-solute  material;  and  this  sorbent, 
in  turn  is,  extracted  using  an  organic  solvent  such  as  hexane.  The  resulting 
concentrated  solute  is  measured  using  a  spect»'ophotometer  and  liquid  chroma¬ 
tograph  .- 

The  column  generator  and  related  high-pressure  liquid  chromatograph  (HPLC) 
methods  of  determining  aqueous  solubility  are  the  most  recent  of  the  three 
described.  However,  It  has  become  the  predominant  method  in  the  literature 
over  the  past  5  years.  In  this  method,  a  column  is  packed  with  an  Inert  solid 
material,  such  as  fine  silica  powder  (100  to  120  mesh),  that  serves  as  a 
support  for  the  organic-solute  phase.  The  fine  powder  provides  a  high  surface 
area  for  the  solute  to  ensure  rapid  equilibration  between  it  and  the  aqueous 
phase.  The  support  is  Impregnated  by  flooding  the  column  with  the  solute  of 
Interest  and  letting  it  sit  for  several  hours.  Water  is  then  admitted  to  push 
the  plug  of  the  organic  phase  from  the  column.  The  column  Is  brought  to  the 
desired  temperature  oy  appropriate  controls  and  flushed  with  water  until  the 
presence  of  undissolved  solute  is  not  detected.  To  ensure  that  equilibrium 
between  the  solute  and  the  aqueous  phase  is  achieved,  the  outlet  concentration 
as  a  function  of  flow  rate  should  be  determined  for  each  compound  of  Interest. 
Ultraviolet,  infrared,  or  spectrofluorescence  detectors  may  be  incorporated 
into  the  HPLC  to  determine  the  solute  concentration  In  each  eluted  phase.  The 
point  at  which  the  solute  concent'^atlons  become  flow-independent  represents 
the  maximum  water  flow  rate  allowed  through  the  column  and  should  also 
represent  the  equilibrium  concentration. 

The  shake-flask  method  has  four  major  sources  of  determining  error 
associated  with  it.  The  first  of  these  stems  from  the  organic-solute  losses 
due  to  adsorption  of  the  solute  on  the  flask  walls  or  to  evaporation.  This 
quantity  Is  somewhat  dependent  on  the  vapor  pressure  of  the  solute,  but  If  the 
quantity  of  the  solution  used  Is  fairly  large  (I.e.,  250  mL),  these  losses  are 
usually  minimized  to  5  percent  of  the  total  solute.  The  second  error  source 
Is  undissolved  solute  suspended  In  the  water  phase.  This  emulsion  may  cause  a 
sample  to  register  an  anomalously  high  solubility.  Centrifugation  or  fil¬ 
tration  should  reduce  this  false  contribution.  The  third  source  arises  from 
the  extraction  efficiency  for  the  removal  of  solute  from  water  using  an  or¬ 
ganic  solvent,  a  technique  common  to  most  methods.  Finally,  the  accuracy  of 
the  detection  method  contributes  to  the  overall  error  In  this  method. 

The  column  method  also  has  several  steps  that  may  result  In  a  negative 
bias  on  the  measured  aqueous  solubility.  The  lack  of  water-solute  equilibrium 
in  the  generator  column  itself  may  produce  an  outlet  concentration  lower  than 
the  "true"  solubility.  Stolzenburg  and  Andren  (Reference  16)  define  this 
certainty  to  approach  95  percent  over  a  tenfold  decrease  In  flow  rates.  The 
efficiencies  for  the  adsorption  and  extraction  steps  and  the  detection 
limitations  may  be  estimated  to  be  similar  to  the  shake-flasx  efficiencies. 


other  problems  are  attendant  to  the  column  method.  If  the  solubility  of 
the  organic  phase  is  greater  than  approximately  1000  ppm,  the  stationary  phase 
on  the  support  depletes  rapidly,  and  a  steady-state  effluent  concentration  may 
not  be  obtained.  A  large  portion  of  chemicals  in  this  study  have  solubilities 
In  and  above  that  range.  Second,  some  strong  UV  adsorbers  (principally 
aromatic  compounds)  may  preferentially  plate  out  on  the  UV  detection  cell 
window.  This  causes  a  significant  signal  drift  during  the  course  of  a  single 
measurement.  Finally,  to  achieve  accurate  calibration  of  the  UV  results  and 
to  avoid  underestimation,  shake-flask  samples  of  known  concentration  must  be 
prepared.-  This  eliminates  most  of  the  proposed  advantages  of  the  generator 
column  technique. 

B.  EXPERIMENTAL 

Initially,  all  three  of  the  principal  measurement  methods  were  examined: 
nephelometry,  generator  column,  and  shake  flask.  The  results  for  some  of  the 
compounds  were  compared  for  the  different  methods.-  By  far,  the  most  satis¬ 
factory  method  turned  out  to  be  the  shake-flask  method.  All  the  final  results 
reported  for  the  solubility  measurements  are  measured  by  this  method.  First, 
the  outcome  of  the  nephelometry  and  generator  column  trials  are  discussed. 

1.  Nephelometry  Measurements 

During  a  preliminary  investigation  designed  to  test  nephelometric 
methods  for  reproducibility  and  consistency  with  literature  values,  it  became 
apparent  that  nephelometry  may  not  be  suitable  for  this  purpose.  Predicted 
values  for  solubility  limits  were  consistently  high  when  compared  to  litera¬ 
ture  values,  experimental  reproducibility  was  poor  at  best,  and  the  relation¬ 
ship  between  light-scattering  and  organic  concentration  appears  to  be  linear 
only  in  a  very  limited  range. 

During  experimentation,  it  appeared  that  loss  of  suspended  organics 
to  the  atmosphere  and  to  glass  surfaces  were  significant  sources  of  error. 
Droplets  of  organic  material  were  observed  to  adhere  to  the  surface  of  the 
glass  sample  containers,  even  while  placed  in  a  sonicating  bath.  Although 
length  of  storage  in  a  closed  container  did  not  (for  nonreactive  components) 
appear  to  alter  eventual  turbidity  readings,  successive  readings  taken  from 
the  same  flask  yielded  decreasing  extents  of  turbidity.  This  indicates  loss 
of  volatile  organics  to  the  atmosphere.  In  addition,  after  short  periods  of 
time  in  a  couvette,  turbidity  gradients  were  evident,  possibly  indicating  the 
tendency  of  dispersed  organic  particles  to  coalesce. 

It  was  also  observed  that  the  assumption  of  detector  linearity  is 
valid  in  a  certain  region  only.  This  is  to  be  expected  because  linearity 
implies  uniform  particle  size,  when  in  reality,  a  size  distribution  exists. 
This  particle-size  distribution  is  a  function  of  the  rate  of  coalescence, 
which  in  turn  depends  on  the  concentration  of  the  dispersed  phase.  In 
addition,  the  distribution  is  dependent  on  the  properties  (e.g.,  surface 
tension)  of  the  organic  phase  in  solution,  so  the  length  of  the  linear  range 
varies  from  one  compound  to  another.  As  a  result,  each  compound  must  be 
examined  independently  before  linearity  can  be  assumed. 
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Preliminary  findings  indicate  that  the  nephelometric  methods  employed 
are  not  sufficiently  accurate  to  warrant  their  use  in  checks  of  experimentally 
determined  solubilities,  in  the  more  dilute  solutions.  However,  the  solubil¬ 
ities  of  several  high  solubility  compounds  on  the  Air  force  list  were  checked 
using  nephelometry.  This  method  was  employed  for  phenol,  toluene,  cis- 
dlchloroethyiene,  trans-dichloroethylene,  methyl  ethyl  ketone,  methyl  isobutyl 
ketone  and  trichloroethylene.  The  solutions  were  prepared  at  several  concen¬ 
trations  well  above  the  anticipated  solubility  limit  in  100  ml  volumetric 
flasks.  These  flasks  were  immersed  in  a  temperature-controlled  ultrasonic 
bath  and  sonicated  until  the  organic  was  completely  dispersed  (2  to  8  hours). 
Samples  were  poured  Into  couvettes,  and  the  turbidity  was  then  determined  by  a 
Hach  2100A  turbidimeter  calibrated  with  standard  (NTU)  solutions.  The  meas¬ 
ured  turbidity  was  plotted  as  a  function  of  the  mixture  concentration  for 
several  dilutions.  The  plot  is  extrapolated  to  zero  turbidity  using  a  least- 
squares  fit,  and  the  concentration  at  that  point  is  used  as  the  solubility 
limit.  This  technique  has  a  fairly  low  accuracy  and  may  only  be  applied 
adequately  to  those  compounds  with  a  solubility  of  over  5000  ppm.  The  nephe¬ 
lometry  results  for  some  of  these  sample  compounds  are  shown  in  Table  13. 
Although  the  method  will  yield  rough  estimates  of  solubilities,  it  is  cer¬ 
tainly  not  accurate  enough  to  predict  changes  in  solubilities  over  small 
(10  ’C)  temperature  intervals. 


TABLE  13.  NEPHELOMETRY  RESULTS 


Compound 

Temperature,  'C 

Solubility, 

ppm 

Correlation 

coefficient 

Toluene 

30 

673.9 

0.989 

20 

601.3 

0.984 

10 

618.9 

0.968 

trans-Oichloroethylene 

30 

6594 

0.970 

20 

6055 

0.998 

10 

5274 

0.995 

cis-Dichloroethylene 

30 

3025 

0.970 

20 

3126 

0.960 

10 

2572 

0.965 

Trichloroethylene 

30 

1282 

0.950 

20 

1250 

0.910 

10 

1271 

0.825 

2.  Generator  Column 

The  major  focus  of  work  performed  in  this  phase  of  the  project  was  on 
determining  the  accuracy  and  reliability  of  the  HPLC-gcnerator  column  method 
for  measuring  solubilities.  Toluene  was  chosen  as  the  Initial  chemical  to  be 
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tested  because  it  has  a  reported  solubility  in  the  median  range  of  the 
chemicals  to  be  analyzed  and  is  a  fairly  common  chemical.  Its  solubility  limit 
has  been  reported  in  several  literature  refe'-ences  and  it  is  also  safe  and 
inexpensive.  The  experimental  procedure  included  testing  various  HPLC  packing 
materials,  varying  flow  rates,  varying  column  loading  techniques,  and  deter¬ 
mining  the  most  accurate  methods  of  calibrating  the  UV  detector.  Testing  was 
done  using  a  Varian  5000  HPLC  unit  equipped  with  a  UV-50  Variable  Wavelength 
Detector. 


The  packing  materials  used  included  treated  and  untreated  Chromosorb 
P  and  glass  beads.-  The  packing  materials  seemed  to  have  little  effect  on  the 
height  of  the  equilibrium  plateau  attained  from  the  analysis.  Changing  the 
flow  rate  from  0.5  cc/minute  to  2.0  cc/minute  a’so  has  no  effect  on  the  signal 
height.  The  equilibrium  value  also  did  not  seem  to  vary  with  the  amount  of 
water  added  to  the  column  when  it  was  being  loaded,  as  long  as  the  concen¬ 
tration  was  well  above  the  solubility  limit.  The  major  problem  encountered 
was  in  calibrating  the  signal  height  to  the  precise  concentration  it  repre¬ 
sents. 


Many  methods  were  used  to  prepare  and  analyze  calibration  standards 
of  known  concentrations.  The  standards  were  prepared  in  vaiying  concentra¬ 
tions  to  obtain  a  linear  response  from  the  UV  detector  of  peak  height  versus 
concentration.  The  results  of  this  study  indicated  that  the  major  drawback  to 
this  method  was  calibration.  The  effort  put  into  producing  reliable  standards 
was  much  greater  than  that  required  to  produce  the  saturateo  solutions.  This 
IS  discussed  in  greater  length  in  the  following  section  for  the  shake-flask 
experiments. 

In  addition,  for  those  compounds  with  solubility  in  excess  of 
approximately  1000  ppm,  a  steady  state  was  difficult  to  achieve  for  long 
periods.  The  UV  detector  response  would  increase  initially  with  the  Intro¬ 
duction  of  the  HPLC  flow  through  the  stationary  phase.  However,  after  a  short 
time  (less  than  5  minutes),  the  UV  response  would  drop  as  the  stationary  phase 
was  depleted.  Larger  columns  saturated  with  the  stationary  phase  might  solve 
this  problem,  however,  the  technique  as  described  in  the  literature  would 
probably  only  be  applicable  to  less  than  half  of  the  test  materials  of  this 
study. 


3.  Shake  Flask  Measurements 

The  primary  method  of  determining  aqueous  solubilities  in  this  study 
was  the  shake-flask  method.  Although  its  concept  is  extremely  simple,  the 
actual  measurements  require  extremely  accurate  laboratory  procedures  to  reduce 
error.  The  shake-flask  method  has  been  applied  to  all  of  the  original 
compounds  except  vinyl  chloride  (gas). 

Calibration  standards  of  known  concentrations  of  the  compound  of 
interest  are  prepared  in  250-milliliter  (+1  mL)  amber  bottles  fitted  with 
double  TeflorP  lined  septa.  The  bottles  are  initially  filled  with  distilled, 
deionized  water  and  sealed  so  as  to  leave  zero  headspace.  The  organic  liquids 
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are  measured  volumetncal ly  and  injected  through  each  septum  using  a  micro- 
liter  syringe.  Conveision  of  parts  per  million  by  volume  to  parts  per  million 
by  weight  is  done  by  multiplying  the  former  by  the  sample  density  (at  25  *C) . 
Solid  samples,  such  as  1 ,4-di chlorobenzene  or  naphthalene,  are  weighed  on  a 
Mettler  microbalance  and  added  to  the  water  before  sealing  the  septum.  The 
saturated  samples  are  prepared  by  adding  organics  far  In  excess  of  the  antici¬ 
pated  solubility  limit  so  that  a  distant  organic  layer  is  formed.  Once  the 
samples  are  prepared,  they  are  shaken  for  1  hour  with  a  wrist-action  shaker 
and  then  allowed  to  rest  until  fully  dissolved  or  equilibrated.  This  equili¬ 
bration  takes  approximately  3  weeks,  on  average. 

The  concentration  measurements  for  the  samples  are  performed  in  two 
ways.  For  those  compounds  that  have  a  strong  UV  adsorbance,  such  as  the 
aromatic  compounds,  the  samples  are  injected  into  a  LOG  UV  III  Monitor  (256 
nm) ,  and  the  UV  adsorbance  for  a  particular  concentration  of  organic  is  deter¬ 
mined.  A  calibration  plot  is  made  of  concentration  versus  adsorbance  for  all 
of  the  independently  prepared  calibration  samples  (4  to  5  samples).  The 
saturated  solution,  which  has  been  equilibrated  at  10  *C,  20  *C,  or  30  *C  is 
then  measured  for  adsorbance  in  the  same  way  as  the  calibration  standards.- 
The  concentration  of  the  saturated  solution  is  then  determined  for  the  extra¬ 
polation  of  the  calibration  line. 

A  second  measurement  method  is  employed  for  those  compounds  that  are 
weak  UV  absorbers.  The  calibration  samples  are  prepared  in  methanol  to 
concentrations  in  the  range  of  the  expected  solubility  limit  in  water.  These 
samples  are  then  injected  directly  into  a  GC  equipped  with  an  80/120  Carbopack 
B/3  percent  SP-1500  column  (1.8  inch  inner  diameter  by  10  feet)  and  a  FID 
detector.  Once  again,  a  calibration  plot  of  GC  response  versus  concentration 
IS  prepared.  A  saturated  aqueous  solution  at  10  ‘C,  20  *C,  or  30  ‘C  is  then 
injected  and  the  GC  response  compared  to  the  plot  to  obtain  the  ultimate 
solubility  concentration. 

The  utility  of  th’S  dual  approach  may  be  demonstrated  for  naphtha¬ 
lene.  Naphthalene  has  a  solubility  in  the  low  ppm  range.  Calibration 
standards  of  1  pom  are  detectable  in  the  gas  chromatograph,  but  extremely 
long  retention  times  (>i  hour)  are  demonstrated  on  any  water-tolerant  columns. 
On  the  other  hand,  narnthalene  has  a  very  strong  UV  adsorbance  in  this  con¬ 
centration  range  and  may  be  measured  accurately  and  rapidly.  Compounds  such  as 
n-hexane  have  virtually  no  UV  adsorption  and  are,  therefore,  not  amenable  to 
this  technique.  Gas  chromatography,  therefore,  becomes  the  method  of  choice. 

The  results  of  the  shake-flask  tests  are  shown  in  Table  14.  This 
table  summarizes  the  aqueous  solubilities  determined  for  the  test  compounds 
for  all  three  temperatures  and  the  method  of  measurement  used  for  each.  A 
detailed  presentation  of  the  calibration  tests  and  95-percent  confidence  bands 
are  shown  in  Appendix  C.  There  were  three  compounds  for  which  solubility 
results  could  not  be  obtained,  most  probably  because  their  solubility  was 
below  the  limits  of  detection:  n-nonane,  decalin,  and  1,2-dinitrotoluene.  In 
addition,  results  for  the  gaseous  compounds,  vinyl  chloride  and  chloroethane 
are  not  available  due  to  the  difficulties  in  assuring  complete  saturation  of 
the  gas  in  water,  as  well  as  the  safety  hazard  associated  in  doing  complete 
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TABLE  14.  LIST  OF  SOLUBILITY  RESULTS 
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saturation  studies.  For  the  sake  of  comparison,  a  compendium  of  litei-ature 
values  for  many  of  the  compounds  is  shown  in  Table  15. 

Comparing  these  literature  values  to  the  ones  obtained  in  the  current 
work,  no  clear  tendencies  (either  high  or  lew)  are  statistically  evident. 
However,  it  should  be  noted  from  the  95  pe'^cent  confidence  bands  listed  in 
Appendix  C  that  In  most  cases  the  temperatures  dependence  is  not  statistically 
significant.  This  is  chiefly  due  to  the  fact  that  solubilities  are  not  a 
strong  function  of  temperature.  Also,  in  examining  the  results  noted  that 
seldom  are  the  solubilities  observed  to  be  monotonic  with  temperature. 

Although  the  effect  of  tsmperatjie  in  these  studies  is  often  not  statistically 
significant,  three  general  trends  may  be  observed.  First,  the  solubility  of 
the  halogenated  hydrocarbons  decreased  with  temperature.  Second,  the  solubil¬ 
ity  of  the  substituted  aromatic  hydrocarbons  increases  with  temperature. 
Finally,  maxima  and  minima  are  observed  for  a  wide  range  of  compounds  without 
any  general  trend  that  can  be  demonstrated  to  be  statistically  significant. 
More  detailed  temperature  studies  for  these  systems  need  to  be  conducted 
beFore  conclusions  as  to  their  origin  may  be  drawn.  A  full  outline  of  the 
statistical  analysis  is  given  n  the  data  analysis  discussion  of  Section  VI. 


TABLE  16.  COMPILATION  OF  AJJUEOUS  SOLUBILITIES  FROM  THE  LITERATURE  FOR  ORGANIC 
COMPOUNDS  OF  INTEREST  TO  THE  AIR  FORCE 


(cont inu*d) 


SECTION  V 


MIXTURE  STUDIES 


In  this  mixture-analysis  phase  of  the  work,  six  components  were  selected 
that,  theoretically,  can  interact  in  multicomponent  solutions:  benzene, 
propylbenzene,  phenol,  1,2-dichlorobenzene,  cis-l,2-dichloroethene,  and  1,2- 
dichloroethane.  These  chemicals  have  structural  features  that  can  cause 
either  self-association  or  the  formation  of  molecular  complexes.  Because 
halogen  substituents  are  electron-withdrawing,  the  chlorinated  compounds 
mentioned  have  induced  dipole  moments  that  are  attracted  to  other  species  with 
similar  part’al  charges.  The  propyl  substituent  in  propylbenzene  also 
exhibits  electron  inductive  effects  and  induced  charges.  With  the  “cis" 
configuration,  the  chlorinated  ethene  compound  has  f-electrons  attracted  by 
the  chlorine  groups  to  produce  a  dipole  moment.  Phenol  is  known  to  self¬ 
associate  by  hydrogen  bonding  in  a  chain  structure,  with  the  length  of  the 
chain  dependent  upon  the  solution  phenol  concentration.  Another  compound  in 
this  study  that  might  interact  with  phenol  is  bis  (2-ethy1hexyl)  phthalate, 
but  detection  limit  problems  encountered  during  the  single-component  runs 
eliminate  it  from  consideration.  Finally,  the  aromatic  compounds  may  "stack 
up"  in  solution  because  of  coordination  of  the  delocalized  »-electrons.- 

The  experimental  matrix  for  the  mixture  runs  was  designed  to  acquire 
qualitative  interaction  data  with  a  minimum  number  of  EPICS  experiments.  This 
is  being  done  with  a  "T"-test  (at  the  95  percent  confidence  level)  of  the 
hypothesis  that  the  Henry's  law  constant  in  the  mixture  is  the  same  as  the 
single-component  value.  Failure  of  the  test  indicates  the  presence  of 
significant  chemical  interactions  in  the  mixture.  The  use  of  several  con¬ 
centration  levels  demonstrates  the  concentration  dependency  of  such  effects, 
if  they  exist.  In  the  mixtures,  the  concentration  of  each  organic  was  diluted 
from  Its  aqueous  saturation  limit  by  a  factor  of  six;  phenol,  however,  was 
also  present  in  one  mixture  at  only  0.5  percent  of  its  saturation  limit  to 
investigate  composition  effects.  Serial  dilutions  of  both  mixtures  are  also 
being  run  as  a  further  independent  check  of  linear  Henry's  law  equilibrium 
behavior. 

The  software  necessary  to  perform  the  statistical  calculations  has  been 
written  in  the  form  of  a  Lotus  1-2-3  spreadsheet/macro  compatible  with  other 
data  analysis  programs  written  for  this  report.  The  T-test  procedure  only 
identifies  significant  differences  in  the  EPICS  average.  However,  if  such 
differences  are  found  to  exist,  the  binary  chemical  interactions  can  be 
quantified  with  an  appropriate  fractional  factorial  experimental  design. 

Given  the  low  concentrations  of  interest  in  this  study,  the  project  team 
originally  believed  that  mixture  effects  would  not  be  observed  and  that  the 
"T"-stat1stic  hypothesis  test  would  confirm  this.  Experimental  work  in  the 
mixture  phase  of  this  investigation,  however,  suggests  the  presence  of 
chemical  interactions  in  solution  even  at  low  dilutions.. 


The  mixture  work  has  focused  on  measurements  at  20  *C  for  five  different 
mixtures  prepared  from  the  six  organic  chemicals  discussed  earlier.  The  test 
samples  were  as  follows: 

1.  Pure  components  in  water 

2.  Mixture  composed  of  equal  volumes  of  saturated  stock  solution  for 
each  organic 

3.  Mixture  identical  to  Mixture  2  except  that  phenol  was  omitted  (a 
volume  of  pure  water  was  added  to  maintain  the  concentrations  of 
the  other  chemicals  at  the  same  level  as  Mixture  2) 

4.,  Mixture  identical  to  HixtU''e  2  except  that  phenol  is  present  at  a 
much  lower  concentration  (e.g.,  a  concentration  of  approximately 
0.5  percent  of  the  saturation  limit) 

5.  Three-component  mixture  consisting  of  benzene,  propylbenzene,  and 
1,2-dichlorobenzene  at  the  same  concentration  levels  as  Mixture  2 

6.  Four-component  mixture  consisting  of  benzene,  propylbenzene,  1,2- 
dichlorobenzene,  and  cis-dichloroethylene  at  the  same  concentra¬ 
tion  levels  as  in  Mixture  2. 

Henry's  law  constants  for  each  component  are  plotted  for  these  mixtures 
(along  with  the  single-component  results  obtained  from  this  study)  in  Figures 
8  through  12.-  The  trends  shown  in  these  plots  are  very  interesting  because 
both  positive  and  negative  deviations  from  ideal  Henry's  law  behavior  are 
indicated.  Comparing  the  results  for  Mixtures  2  and  4,  an  increasing  phenol 
concentration  in  multicomponent  solution  raises  the  air-water  partition  co¬ 
efficient  for  four  of  the  other  components:  benzene,  propylbenzene,  cis- 
dichloroethylene,  and  1,2-dichloroethane.  Apparently,  strong  phenol  self¬ 
association  at  higher  concentration  reduces  the  affinity  of  these  compounds 
for  the  liquid  phase.  The  final  chemical,  l,2-dichlorobenz8ne,  did  not  appear 
to  be  affected  significantly  by  changes  in  the  phenol  concentration  and  was 
therefore,  included  in  all  five  mixtures  as  a  "benchmark"  compound  to  confirm 
that  the  EPICS  technique  was  giving  accurate  results. 

An  observed  partition  coefficient  that  is  lower  than  the  true  Henry's  law 
constant  represents  a  negative  deviation  from  ideal  behavior.  For  benzene  and 
propylbenzene,  however,  positive  deviations  from  ideal  solution  behavior  were 
also  noted  during  the  mixture  tests.  The  partition  coefficients  measured  for 
these  two  chemicals  in  Mixture  3  (no  phenol)  were  significantly  higher  than 
the  Henry's  law  constants  determined  in  the  single-component  measurements,  in¬ 
dicating  a  mixture  effect  not  attributable  to  the  phenol  concentration.  It 
appears  to  be  a  “saltiiig-out"  effect  in  which  one  or  more  of  the  other  or¬ 
ganics  present  in  the  mixture  increase  the  hydrophobicity  of  benzene  and  pro¬ 
pylbenzene  (e.g.,  lowers  their  respective  aqueous  concentration). 
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Mixture  Study  Results  for  Benzene. 
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Figure  10.  Mixture  Study  Results  for  1,2-Dichloroethane. 


Figure  11.  Mixture  Study  Results  ror  Cis-Dichloroethylene. 
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Mixture  Study  Results  for  1,2'Dichlorobenzene. 


After  observing  the  above  trends,  serial  dilutions  Mixtures  2,  3,  and  4 
were  analyzed  to  verify  any  deviations  from  ideal  Henry's  law  behavior.  Sur¬ 
prisingly,  all  compounds  except  phenol  were  found  to  have  a  linear  relation¬ 
ship  of  vapor  composition  to  liquid  composition  over  the  entire  dilution 
range.  Linear  regressions  of  the  mixture  serial  dilution  data  resulted  in 
correlation  coefficients  (r  squared  values)  better  than  0.99  for  the  mixture 
compounds.  Despite  the  linearity  of  the  dilution  curves,  the  project  team 
believes  that  the  deviations  from  ideal  behavior  indicated  by  the  EPICS  tech¬ 
nique  are  real  and  reproducible. 
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SECTION  VI 


MODELING 


A.  COMPUTER  DATA  ACQUISITION  SYSTEM 

Extensive  computer  software  was  developed  to  satisfy  the  data  manipulation 
requirements  of  the  current  Henry's  law  constant  measurement  program.  The 
EPICS  phase  of  the  effort,  in  particular,  has  involved  extensive  "bookkeeping" 
and  data  reduction.  To  address  this  need,  a  software  package  composed  of 
BASIC  and  LOTUS  1-2-3  programs  was  written  to  completely  automate  the  data 
acquisition  process. 

The  first  step  in  the  EPICS  data  collection  procedure  is  to  transmit  the 
integrated  peak  a-ea  from  a  Hewlett-Packard  (HP)  3392A  integrator  to  a  per¬ 
sonal  computer  (an  IBM-compatible  Tandy  1200  HD  system)  in  ASCII  format  by  way 
of  an  RS-232C  interface.  It  has  also  proven  desirable  to  provide  the  inte¬ 
grator  with  bookkeeping  iiiforraation  from  the  computer  keyboard.-  This  is  pos¬ 
sible  because  the  HP  3392f,  with  its  own  language  syntax,  can  conduct  two-way 
"conversations"  with  computers  by  sending  and  receiving  appropriate  ASCII 
character  strings.  A  BASIC  program  written  to  operate  tie  asynchronous  com¬ 
munication  link  between  the  two  devices  makes  use  of  the  HP  3392A  language 
convention  for  information  exchange. 

Using  this  capability,  the  laboratory  analyst  designates  information,  such 
as  a  component  identifier,  the  run  (bottle)  equilibration  temperature,  the 
replicate  number,  etc.,  that  will  be  inserted  onto  the  printer/plotter  copy  of 
the  report  and  permanently  stored  in  a  computer  report  file.  The  descriptive 
run  information,  when  cross-referenced  against  a  properly  kept  laboratory 
notebook,  is  an  excellent  precaution  that  helps  ensure  data  integrity.  Termi¬ 
nation  of  a  given  run  by  the  analyst  at  the  integrator  console  causes  immedi¬ 
ate  transmission  of  the  peak  area  report  to  the  computer.  In  addition,  if  a 
previous  report  is  stored  in  the  integrator  memory  (active  workspace),  it  may 
be  "customized"  and  retransmitted. 

After  the  report  is  stored  in  a  "capture"  file  by  the  BASIC  program,  LOTUS 
1-2-3  routines  developed  specifically  for  this  application  are  then  used  to 
scrutinize  the  c  ,a.  The  testing  order  was  randomized  for  statistical  valid¬ 
ity,  but  the  data  manipulation  was  made  more  complex.  The  LOTUS  routines  are 
therefore  designed  to  rearrange  the  data  into  a  logical  organizational  pattern 
with  descriptive,  unique  names  for  the  data  files. 

For  instance,  one  of  the  LOTUS  programs,  a  "keyboard  macro,"  (in  LOTUS 
terminology)  is  used  to  prccass  and  organize  an  entire  day's  integrator  re¬ 
ports.  The  routine  does  this  by  examining  a  capture  file  line  by  line,  brac¬ 
keting  each  report  after  detecting  the  beginning  run  number  and  the  end-of- 
report  character,  and  then  extracting  the  pertinent  peak  areas  by  position.  A 


typical  capture  file  contains  all  runs  for  1  flay  of  analysis,  with  each  run 
consisting  of  four  injections  for  a  single  rompound  at  a  specified  tempera¬ 
ture.  Each  set  of  extracted  component  data  is  given  a  file  name  derived  from 
the  component  identification  number,  the  temperature  code,  and  the  replicate 
number. 

Another  LOTUS  spreadsheet  then  accesses  these  individual  data  files  on 
request  and  performs  the  necessary  calculations  and  statistical  tests.  One 
such  analysis  program  is  created  for  each  component  and  can  be  updated  as  new 
information  is  received.  Calculated  statistical  Quantities  include  the  co¬ 
efficient  of  variation  relative  standard  deviation)  for  all  replicate  Henry's 
law  constant  observations,  the  leraperature  regression  parameters  (slope  and  y- 
intercept)  and  associated  correlation  coefficient  for  each  component,  and  the 
Student's  "T“-test  confidence  bands  for  both  the  raw  data  and  the  temperature 
regression  predictions.  Using  these  statistics,  it  was  possible  to  determine 
the  overall  precision  of  the  EPICS  procedure  for  all  components  and  to  esti¬ 
mate  the  accuracy  bounds.  Appropriate  temperature  regression  and  confidence 
interval  plots  were  generated  automatically  by  the  software  to  present  trends 
graphically.  A  final  LOTUS  program  was  written  to  combine  the  data  for  all  oi 
the  chemicals  into  a  master  tabulated  summary  convenient  to  view  or  mani¬ 
pulate. 

Divided  according  to  component,  the  information  generated  by  the  LOTUS 
software  for  each  chemical  is  composed  of  the  following:' 

•  Two-page  tabulation  of  the  injection  peak  areas,  Henry's  law 
constant  estimates,  and  Coefficient  of  Variation  COV  for  the 
component  at  give  temperatures  (10,  15,  20,  25,  and  30  *C) 

•  Temperature  regression  plot  (In  H  versus  1/T) 

•  Plot  showing  the  95  percent  confident  band  on  the  temperature 
regression  predictions 

•  Plot  illustrating  the  95  percent  limits  (lower  and  upper)  on  the 
averages  of  the  estimates  calculated  at  each  temperature. 

All  of  this  information  is  displayed  for  the  48  compounds  of  interest  in 
Appendix  8  of  this  report. 

Similar  computer  programs  have  been  developed  for  the  bubble  column  and 
solubility  tests.  Aga.n,  the  emphasis  was  on  ease  of  data  reduction  and  mani¬ 
pulation,  plus  valid  statistical  analyses.  In  summary,  all  software  items 
used  in  the  current  program  allow: 

•  User-friendly  operation 

•  Simplified  data  handling/reduction 

•  Valid  statistical  tests  and 

•  Elimination  of  human  error. 
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Because  the  program  functions  are  completely  menu-driven  rapid  data  processing 
by  the  laboratory  analyst  is  possible,  thereby,  providing  useful  feedback  for 
technique  improvement.  Most  important  is  the  removal  of  the  human  element 
during  the  data  analysis,  which  has  dual  benefits  of  eliminating  manual  log¬ 
ging  errors  and  ensuring  computational  accuracy., 

B.  CORRELATION  DEVELOPMENT 

To  predict  the  thennophysical  properties  of  VOCs  in  water  a  thermodynamic 
framework  for  the  behavior  of  pure  compounds  in  water  is  needed.  Three  dif¬ 
ferent  thermodynamic  techniques  for  correlating  experimental  Henry's  law  con¬ 
stants  and  aqueous  solubilities  have  been  examined.  Two  of  these  techniques 
are  based  on  modeling  the  VOCs  as  collections  of  discrete  molecular  fragments: 
the  UNIQUAC  Functional  group  Activity  Coefficient  (UNIFAC)  model  of  Fredens- 
lund  et  al.  {Reference  3)  and  the  substituent  constant  method  introduced  by 
Fujita  et  al.  (Reference  42)  and  described  in  detail  by  Leo  et  al.  (Reference 
4).  The  third  technique  is  based  on  a  combined  chemical  and  physical  theory 
of  thermodynamic  properties,  the  so-called  associated  solution  theory,  as 
recently  applied  to  Henry's  law  constant  of  supercritical  gases  by  Hu  et  al . 
(Reference  43).  The  purpose  of  using  these  multiple  techniques  is  to  cri¬ 
tically  evaluate  three  different,  but  valid,  approaches  for  their  appli¬ 
cability  to  environmental  systems  and  their  predictive  capacity  for  unmeasured 
multicomponent  systems. 

1.  Unifac  Methods 

a.  Basis  of  UNIFAC 

The  UNIFAC  model  was  developed  in  1975  to  correlate  large  quan¬ 
tities  of  data  with  a  few  molecular  parameters.  The  first  group  contribution 
technique  developed  for  activity  coefficient  prediction  was  that  of  Wilson  and 
Deal  (Reference  12);  the  subsequent  development  of  UNIFAC  by  Fredenslund  et 
al.  (Reference  3)  owes  much  to  their  early  work.  Activity  coefficients  are 
closely  related  to  various  excess  properties.  In  particular,  the  excess 
G’bbs'  free  energy,  G-^',  may  be  related  to  the  activity  coefficient  in  a 
straight  forward  manner: 


6!"=  RT  In  7*  .  (21) 

In  the  development  of  the  UNIQUAC  equation  for  phase  equilibria  Abrams  and 
Prausnitz  (Reference  44)  used  the  the  concept  of  local  chemical  compositions 
to  stipulate  that  the  excess  Gibbs'  free  energy  is  made  up  of  two  parts:' 


•  A  contribution  due  to  the  difference  in  volume  and  surface  area 
of  molecules  (a  configurational  part) 

•  A  contribution  due  to  the  energetic  interactions  between  molecules  (a 
residual  part). 


In  terms  of  the  activity  coefficient,  this  may  be  written: 
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In7,  =  ln7,co''f'g“''ational  +  107, residual  .  (22) 

The  parameters  that  appear  in  UNIFAC  are  R|<  and  Q)(,  which  are  van  der  Waals 
group  volumes  and  group  surface  areas,  respectively,  a  molecular  coordination 
number,  z;  and  two  binary  local  composition  group  interaction  parameters,  A,j 
ana  Aj,.  The  coordination  number  is  arbitrarily  set  to  10,  and  then  the  van 
der  Waals  group  volumes  and  group  surface  areas  are  obtained  directly  by 
rescaling  Bondi's  (Reference  45)  tabulated  values.  The  binary  interaction 
parameters  have  been  obtained  for  more  than  40  different  groups  by  fitting  the 
model  to  some  55,000  data  points,  including  vapor-liquid  equilibria  (VLE) , 
liquid-liquid  equilibria,  infinite  dilution  activity  coefficients  and  azeo¬ 
tropic  points  (Reference  13).  Of  greatest  interest  to  environmentalists  may 
be  the  recent  compilation  of  infinite  dilution  activity  coefficients  (Refer¬ 
ence  46)  and  its  application  to  upgrading  the  UNIFAC  parameters  for  vapor- 
liquid  equilibria. 

Several  variations  of  UNIFAC  have  been  developed  recently.  The  dif¬ 
ferences  in  the  various  models  include  changes  in  the  combinatorial  contri¬ 
bution  to  the  excess  free  energy,  alternate  choices  for  subdividing  molecules 
into  groups,  the  introduction  of  temperature-dependent  parameters,  and  inclu¬ 
sion  of  solvation  and  chemical  association  effects.  The  major  difficulty  in 
implementing  these  variants  on  UNIFAC  is  the  lack  of  a  fitted  parameter  base. 
The  original  UNIFAC  model  with  the  latest  VLE  parameter  base  was  used  as  the 
initial  benchmark  calculation  against  which  to  compare  our  experimental 
Henry's  law  constants.  The  combinatorial  contribution  to  UNIFAC  can  be 
modified,  as  suggested  by  Kikic  et  al.  (Reference  47).  The  binary  UNIFAC 
parameters  can  then  be  refitted  in  this  same  manner  to  our  data  base  of  ex¬ 
perimentally  obtained  infinite  dilution  activity  coefficients.  In  addition, 
these  group  interaction  parameters  may  be  fittea  to  liquid-liquid  solubilities 
of  VOCs  in  water..  Then  the  extent  to  which  optimized  parameter  values  differ 
between  the  two  data  sets  can  be  determined.  This  modification  of  UNIFAC 
offers  a  potentially  more  accurate  method  of  predicting  Henry’s  law  constants 
for  environmental  systems  than  the  original  version  of  UNIFAC. 

b.  UNIFAC  Parameter  Fitting 

To  satisfy  tne  projects  thermodynamic  modeling  objectives,  UNIFAC  bi¬ 
nary  interaction  parameters  were  fitted  to  the  experimental  VLE  (Henry's  law 
constant)  data  for  the  aliphatic  and  aromatic  chemicals  using  a  simplex  algo¬ 
rithm  written  for  this  purpose.  A  sequential  fitting  procedure  was  used 
rather  than  simultaneously  fitting  all  the  interaction  parameters  covered  by 
the  experlmen'a I  data.  The  sequential  method  required  much  less  conputer 
time,  and  was  judged  by  the  project  team  to  be  better  suited  for  updating  the 
correlation  database  at  this  time;  however,  a  simultaneous  fitting  algorithm 
may  be  used  in  the  future. 

Ihe  basic  objective  was  to  model  the  liquid  phase  activity  coeffi¬ 
cients  of  volatile  organics  infinitely  dilute  in  water  using  UNIFAC.  Then, 
given  experimental  Henry's  law  constants,  infinite  dilution  activity  coef¬ 
ficients  and  optimal  binary  Interaction  parameters  for  UNIFAC  were  calculated. 


A  program  was  written  to  calculate  infinite  dilution  activity  co¬ 
efficients  from  experimentally  determined  Henry's  law  constants  (H,)  for  the 
aliphatic  and  aromatic  chemicals  involved  in  this  study.  The  temperature 
regression  and  coefficient  of  variation  results  for  these  compounds  were  very 
good,  with  only  a  few  data  points  omitted  after  a  visual  examination  of  the 
data  quality.  To  perform  these  calculations,  pure  component  vapor  pressure 
data  and  density  data  for  water  at  the  temperatures  of  interest  were  required. 
A  search  was  conducted  for  vapor  pressure  constants  for  the  Antoine  Equation 
and  the  Harlacher-Braum  (HB)  equation.  Harlacher-Braum  constants  are  scarce, 
and  solution  of  the  H-B  equation  requires  implementation  of  a  Newton-Raphson 
iteration  scheme,  initial  guesses,  etc.  Because  the  Antoine  equation  predic¬ 
tions  are  within  a  few  percent  of  the  H-B  predictions,  the  simpler  Antoine 
equation  was  implemented. 

The  fugacity  coefficient,  involved  in  the  calculations  of  7*00  ^1  ■ 

was  plotted  against  absolute  temperature,  T.  Over  the  temperature  range  10- 
30  *C,  this  quantity  was  wel 1 -approximated  by  a  line,  enabling  simple  linear 
interpolation  of  discrete  data  (see  Figure  13).  A  copy  of  the  FORTRAN  code, 
as  well  as  sample  input  and  output  files  to  perform  the  calculation  of  Tvqc 
from  Hi  may  be  found  in  Appendix  E  (Volume  II).  The  input  file  SUHHARY.DAT  is 
a  file  of  the  experimental  Henry's  law  constants,  temperatures,  etc.  the 
numbering  of  which  does  not  exactly  correspond  to  the  numbering  of  the  vapor 
pressure  data  set  (ANT0INE.DAT).  The  compound  list  used  in  the  vapor  pressure 
data  set  is  also  found  in  Appendix  E.  Only  the  aliphatic  and  aromatic 
chemicals  in  the  present  work  were  used  in  the  development  of  the 
environmental  UNIFAC  correlation. 

(1)  Aliphatics 

To  test  the  approach  for  the  data  analysis,  aliphatic  water  binary 
interaction  parameters  were  determined,  based  on  data  for  4  compounds:  n- 
nonane,  n-hexane,  2-methylpentane,  and  cyclohexane.  Six  different  objective 
functions  were  tested  with  the  optimization  routine  to  determine  the  best 
functional  form  to  minimize  when  fitting  interaction  parameters.  Tnese 
results  are  presented  in  Table  16.  The  objective  function  chosen  was  the 
relative  difference  in  the  natural  logarithms  of  the  experimental  and  cal¬ 
culated  activity  coefficients.  This  has  a  sound  theoretical  basis,  as  may  be 
seen  from  the  Gibbs'  free  energy  relationship  In  Equation  (21).  The  experi¬ 
mental  data  on  which  this  fit  is  based,  and  the  comparison  of  experimental, 
and  predicted  'nfinite  dilution  activity  coefficients  (as  well  as  those 
predicted  by  the  existing  VLE  database)  are  shown  in  Figures  13  through  16. 
These  plots  show  that  the  new  L'HIFAC  data  base  provides  a  very  large  improve¬ 
ment  in  activity  coefficient  predictions.  The  improvement  demonstrates  the 
merit  of  this  VLE  modelling  approach,. 


(2)  Aromatics 

To  optimize  the  binary  interaction  parameters  associated  with 
aromatic  compounds,  a  much  more  involved  procedure  is  required.  To  Improve 
the  parameter  fit,  experimental  activity  coefficients  were  plotted  against 
temperature  for  each  compound  and  obvious  outlying  points  were  excluded  after 
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Figure  13.  Infinite  Dilution  Activity  Coefficients  for  2-Methyl  Pentane. 

!  I  -  Experimental,  +  -  Old  VLE  Binary  Interaction  Parameters, 
O  -  New  Parameters  from  This  Study. 
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Figure  15.  infinite  Dilution  Activity  Coefficients  for  n-Nonane. 

i  I  -  Experimental,  +  -  Old  VLE  Binary  Interaction  Parameters, 
O  -  New  Parameters  from  This  Study. 
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Figure  16.  Infinite  Dilution  Activity  Coefficients  for  n-Hexane. 

I  I  -  Experimental,  +  -  Old  VLE  Binary  Interaction  Parameters 
O  -  New  Parameters  from  This  Study. 


TABLE  16.  LIST  OF  OBJECTIVE  TRIAL  FUNCTIONS  FOR  UNIFAC 
PARAMETER  OPTIMIZATION 


Objective  function 


Binary  interaction  parameters 


Ah20-CH3 


AcH3-H20 


Original  VLE 


''^'calc  -  T^exptl 
X,  =  0 


300 

738 


1318 

0.0 


^^calc  *  ^exptl^ 
X,  =  1  ppm 


738 


0.07 


Tcalc  ■  Wl) 


X^  =  1  ppm 


484 


5280 


“^calc*  '^exptl 
’'^exptl 


745 


2.6 


-ycaic-  WlV 


498 


3523 


'ycalc-  Wl' 

.  ’^'exptl 

'I'exptl  =  0.99999 


493 


4168 


visual  inspection.  The  previously  optimized  Ah^o-CH3  interaction  parameter 
was  included  in  ihe  VLE  data  set.  Because  7“  is  insensitive  to  changes  in 
Ach3-H20i  **45  not  modified. 

The  interaction  parameter  Ah20-ACH  **4S  determined  based  on  the 
data  for  benzene  only.  This  value  was  substituted  into  the  VLE  data  set,  and 
the  new  interaction  parameters  -  Aujo-accHi  Aach-acch3i  snd  Aacch3-ach  ~  were 
then  determined  based  on  data  for  orcho,  meta,  and  para-xylene;  trimethyl- 
benzene;  and  toluene.  The  VLE  data  set  was  again  updated  to  reflect  changes 
in  these  parameters.  To  continue  the  sequence,  a  new  value  of  An20-ACH  was 
determined  based  on  the  data  foi  benzene,  the  xylenes,  trimethyl  benzene,  and 
toluene,,  and  database  again  updated.  Finally,  new  values  for  the  interaction 


parameters  -  Ah20-ACi.H3.  AaCH-ACCH3.  and  AAccH3-ACH-were  determined  from  the 
revised  data. 

This  sequence  of  updating  ana  refitting  continues  iteratively 
until  the  the  numerical  values  of  each  of  the  five  interaction  parameters 
converges  to  a  constant  value.  From  trial  runs,  this  solution  does  indeed 
appear  to-  be  converging  with  each  iteration.  Results  from  this  procedure  are 
shown  In  Table  17.  The  next  portion  of  the  overall  modeling  effort  would  be 
to  increase  the  scope  of  this  iteration  series  to  Include  alkyl-substituted 
aromatics  {e.g.,  methyl  ethylbenzene),  and,  finally,  the  halogenated  com¬ 
pounds.  Unfortunately,  this  expanded  fitting  scheme  was  beyond  the  scope  of 
the  current  work. 

TABLE  17.  FITTING  PROCEDURE  FOR  BINARY  INTERACTION  PARAMETERS 


*H20-CH3  =  492.6  (from  optimization  of  aliphatic  hydrocarbon/water  data) 


Original  VLF 
data  base 

Iteration 

1 

Iteration 

2 

AH20-ACH  =  362.3 

533.218 

(based  on  benzene) 

443.966 

AH20-ACCH3  =  377.6 

133.554 

214.480 

%H-ACCH3  =  167  0 

-300.469 

-604.445 

AACCH3-ACH  =  -146  8 

278.024 

447.705 

2.  Group  Contribution  Metnod 

Tne  second  thermodynamic  technique  investigated  for  correlating 
Henry's  law  constants  of  VOCs  in  water  is  the  simple  substituent  constant 
method  of  Leo  et  a1.  (Reference  4).  By  analogy  with  the  Hammett  constant 
technique  from  physical  organic  chemistry,  Fujita,  et  al.  (Reference  42) 
developed  a  substituent  constant  approach  for  computing  Henry's  law  constants 
and  octanol-water  partition  coefficients,  which  relies  on  the  additivity  of 
the  free  energetic  contributions  to  both  properties.  This  technique  is 
described  in  full  in  Leo  et  al.  (Reference  4)  and  has  some  potential  advan¬ 
tages  over  tne  UNIFAC  method.  Despite  its  strictly  correlational  approach,.  It 
can  in  principal  account  for  inductive,  resonant,  sterlc,  branching,  and 
conformational  effects  that  are  completely  ignored  by  the  UNIFAC  model. 

The  correlation  of  equilibrium  air-'water  partition  coefficients  with 
molecular  structure  is  a  fonnidable  task  due  to  the  difficulty  in  isolating 
the  competing  functional  group  influences.  Attempts  to  predict  partition 
coefficients  using  a  group  contribution  approach  have  proven  successful  when 
the  reference  system(s)  have  been  chosen  so  that  similar  steric  and  electronic 


interactions  are  present  in  the  derivative  compounds.  The  additivity  constant 
(t)  for  a  given  substituent  group  (designated  oy  "X")  is  defined  by  Leo  et  al. 
(Reference  4)  as  follows: 


(23) 


where 


Px  =  partition  coefficient  of  the  derivative  compound 

=  partition  coefficient  of  the  parent  molecule 

=  logarithm  of  the  partition  coefficient  for  the  function  X. 

Summing  the  appropriate  Tx  values  for  the  various  functional  group  and  struc¬ 
tural  elements  of  a  compound  of  interest  gives  the  logarithm  of  its  predicted 
partition  coefficient. 


Values  of  r  correlated  for  a  substituted  species  in  a  homologous 
series  have  been  shown  by  Leo  et  al.  (Reference  4)  to  be  relatively  constant 
for  the  octanol-water  system,  for  example,  the  addition  of  methyl  groups  to 
various  benrene  derivatives  has  been  found  to  yield  t-ch3  values,  which  are 
surprisingly  constant  at  around  0.5. 


There  is  thermodynamic  justification  for  assuming  that  Henry's  law 
constant  may  obey  the  form  of  Equation  (23),  if  one  considers  Henry's  law 
constant  to  be  analogous  to  an  equilibrium  constant.  Letting  N  denote  th  ■ 
parent  molecule  and  X  the  fragment  of  interest,  the  following  reactions 
describe  the  phase  equilibria  between  air  and  water  for  the  parent  and  its 
derivative: 


PARENT 


(N-H) 


H.O 


(N-H) 


air,  ■'1 


TliTH)^] 


K. 

(N-x)  ; 

DERIVATIVE  H^O 

where 


TTiTxyi^] 


(25) 


Ki  =  equilibrium  constant 
(]  =  species  concentration. 


The  quantity  from  the  form  of  Equation  (23),  is  therefore 


T 


logjo 


(26) 
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with  (K2/K1)  being  the  equi  1  ibriiim  constant  for  the  overall  reaction:’ 


«^-«)air‘  <N-X)'a,,MN-H^o  .  (27), 

The  free  energy  change  resulting  from  the  introduction  of  the  substituent  X  is 
obtained  from  this  reaction  as  follows: 

AGtotal  =  ’°9l0  -  «{N-X)  ^  «(N-H)  -  <'{N-X)  ^ 

The  assumption  that  the  free  energy  of  an  individual  molecule  can  be 
represented  as  the  sum  of  the  free  energies  of  its  part  plus  their  interac¬ 
tions  gives  rise  to  the  expressions:- 


‘^(N-H)  =  S  S  Sh 


‘^(N-X)  =  S  ^  ®X  *  Sx  • 

Putting  these  equations  into  Equation  (15)  gives  the  result: 
2.303  RT  logjg  ;(K2/Kj)  =  AG^  +  -  AG^j^  -  AG„ 


„air  ~H_0 

AGj  =  Gj  ■ 


If  the  interaction  terms  may  be  neglected,  we  get  the  following  ex¬ 
pression  for  ir^: 

2.303  RT  logjg  (K^/Kj)  =  AGj^  -  AG^  (32) 


AGy  -  AGm 

*^x  "  "  2.303  RT 


The  definition  of  as  shown  here,  is  physically  defensible  because  it  is 
related  to  the  free  energy  change  in  a  partitioned  compound  that  results  when 
a  substituent  X  replaces  a  hydrogen  atom  of  the  parent  compound. 

3.  Intramolecular  Influences  on  Partitioning 

For  development  of  a  comprehensive  fragment  correlation,  it  is  neces¬ 
sary  to  isolate  and  quantify  at  least  five  different  types  of  functional  group 
effects  on  Henry's  law  constant: 


Inductive  effects 


•  Effects  of  resonance 

•  Steric  effects 

•  Effects  of  branching 

•  Conformational  effects  and  "rigidity." 

The  following  paragraphs  discuss  these  structural  Influences  and  give 
qualitative  examples  of  their  relatlonshlo  to  partitioning  behavior. 

a.  Inductive  Effects 

Inductive  effects  are  the  result  of  electron-withdrawing  sub¬ 
stituents  that  make  the  lone-pair  electrons  belonging  to  a  second  functional 
group  fsuch  as  OH)  less  available  for  hydrogen  bonding.  This  decreases  the 
affinity  of  the  compound  tor  the  water  phase,  thus,  raising  the  value  of 
Henry's  law  constant.  A  group  contribution  correlation  can  account  for  this 
behavior  by  including  x  values  for  the  addition  of  electron-withdrawing  groups 
(such  as  NO2)  to  parent  compounds  containing  functional  groups  with  electrons 
available  for  hydrogen  bonding. 

b.  Effects  of  Resonance 

The  effect  of  electron  delocalization  on  r  values  should  also  be 
included  in  a  group  contribution  correlation.  For  octanol-water  systems,  the 
partition  coefficient  Increases  as  the  result  of  iransferring  a  functional 
group  from  an  aliphatic  to  an  aromatic  position.  Similar  behavior  may  also  be 
anticipated  for  Henry's  law  constant.  Another  electronic  effect  worth  study¬ 
ing  IS  the  presence  of  multiple  bonds  (double,  triple)  as  well  as  conjugated 
double  bonds. 

c.  Steric  Effects 

Steric  effects  are  quite  varied  In  nature,  with  the  most  common 
perhaps  being  a  case  in  which  the  Henry's  law  constant  Increases  because  lone- 
pair  electrons  (on  an  OH  group,  for  nstance)  are  shielded  from  the  sur¬ 
rounding  water  molecules.  This  blocking  of  a  hydrophyllc  functional  group 
could  be  the  result  of  Inert  alkyl  groups,  crowding  of  functional  groups,  or 
chain  branching.  Intuitively,  one  would  expect  that  all  such  circumstances 
would  decrease  the  hydrogen  bonding  capability  of  the  compound. 

d.  Effects  of  Branching 

Branching  Is  another  Important  factor  In  correlating  Henry's  law 
constant  with  chemical  structure.  Work  with  the  octanol-water  system  has 
already  shown  (Reference  4)  that  branching  results  in  lower  substituent  x 
values  than  for  a  "normal "  configuration  when  aliphatic  compounds  are  con¬ 
sidered.  Furthermore,  when  branching  in  that  chemical  class  occurs  at  the 


functional  group,  the  effect  appears  to  be  slightly  greeter,  with  the  contrast 
between  the  normal  and  "branched”  configurations  becoming  more  evident. 

Because  of  uncertainties  in  this  area,  however,  a  constant  value  of  r  for  the 
effect  of  branching  on  the  octanol-water  partition  coefficient  has  been  chosen 
(Reference  4),  and  it  may  prove  desirable  to  take  the  same  approach  for 
Henry's  law  constant. 

e.  Conformational  Effects 

Conformational  effects  are  closely  related  to  rigidity  and  both 
should  be  considered  in  a  structural  correlation.  For  example,  a  simple  case 
of  conformational  influences  occurs  when  long  aliphatic  chains  tend  to  coil  up 
on  solution,  forming  molecular  oil  droplets.  Another  very  interesting  example 
occurs  in  certain  phenyl  derivatives  when  the  dipole  of  a  side  chain  interacts 
with  the  delocalized  electrons  of  the  ring.  This  .auses  "folding"  of  the  side 
chain  onto  the  phenyl  ring  that  promotes  Intramolecular  hydrogen  bonding. 
Presumably,  such  attractions  would  tend  to  give  vhe  compound  an  aqueous  solu¬ 
bility  greater  than  expected,  thus,  towering  the  anticipated  Henry's  law  con¬ 
stant. 

f.  Rigidity  Effects 

The  effect  of  the-  rigidity  of  a  molecule  on  tne  logarithn  of  its 
Henry's  law  constant  (a  quantity  analogous  to  1  ipophi 1 icity)  may  also  be 
important.  Because  the  Henry's  law  constant  of  a  compound  will  be  calculated 
from  small,  flexible  fragments,  restricted  rotation  In  the  compound  may  have 
an  adverse  effect  on  the  accuracy  of  the  correlation.  Based  on  the  ideas  of 
Aranow  and  Witten  (Reference  48),  rigid  molecules  should  have  a  lower  experi¬ 
mental  Henry's  law  constant  than  that  predicted  by  such  a  correlation.  Ihis 
suggests  two  conclusions:- 

♦  A  ''rigidity  factor"  should  be  Included  in  a  preaicti»e  corre¬ 
lation. 

•  Such  a  factor  would  be  entropy  related  and,  thus,  temperature 
dependent. 

The  second  conclusion  'mplies  that  If  the  rigidity  changes  significantly  in  a 
homologous  series  of  compounds,  the  effect  on  the  temperature  dependency  of 
the  structural  contributions  to  Henry's  law  constant  could  be  significant. 

4.  Other  Approaches 

The  final  thermodynamic  technique  evaluated  1r  tills  study  is  the  as¬ 
sociated  solution  theory.  Originally  developed  as  a  purely  chemical  theory  of 
nonidealities  by  Oolezelak  in  1908,  (Reference  49).  the  method  has  been  re¬ 
cently  compined  with  physical  theories  of  jnldealit’"-  to  generate  models  of 
Henry's  law  constants  in  water  (Reference  43).  Thy  v.  !e1  Is  based  on  a  four- 
step  calculation  of  the  free  energy  of  dissolving  a  gaseous  component  in 
water:  (;)  formation  of  pure  associated  and  solvated  species  in  the  perfect 
gas  stale,  (2)  mixing  of  the  pure  associated  and  solvated  species  in  the 


perfect  gas  state,  (3)  introducing  a  hard-sphere  contribution  to  the  energy, 
and  (4)  imposing  an  unattractive  potential  on  the  hard-sphere  mixture  of 
associated  and  solvated  species.  The  model  has  proven  successful  in  corre¬ 
lating  Henry's  law  constants  of  pentane,  hexane,  and  benzene  in  water  at 
temperatures  from  0  to  300  *C.  Future  correlation  work  may  include  testing 
the  model  against  the  existing  environmental  data  oase  and  extending  it  to 
mixtures  for  comparison  with  the  modified  UNIFAC  method  and  the  Hansch  sub¬ 
stituent  group  method. 

C.  COMPUTER  SOFTWARE 

The  UNIFAC  algorithm  has  been  implemented  as  a  FORTRAN  computer  program 
that  can  be  run  on  any  IBM-compatible  (MS-DOS)  personal  computer.  'With  this 
software,  Henry's  law  constants  and  aqueous  solubilities  can  be  calculated  for 
"new"  compounds  (i.e.,  those  not  represented  in  the  current  VLE  database)  from 
their  constituent  structural  grouos.  Partition  coefficients  for  organic 
compounds  in  multicomponent  mixtures  can  also  be  calculated  with  the  UNIFAC 
software.  The  simplex  algorithm  for  fitting  binary  interaction  parameters  to 
experimental  data  in  the  manner  explained  above  is  included  in  the  software 
package  to  give  the  Air  Force  the  ability  to  update  the  UNIFAC  database  with 
future  in-house  measurements.  At  present,  this  computer  simplex  routine  is  in 
"batch"  form  (i.e.,  input/output  operations  are  not  interactive). 

To  make  the  VLE  and  LLE  algorithms  easy  to  use,  the  FORTRAN  routines  are 
written  in  modular  form,  and  the  modules  connected  by  interactive  menus  and 
graphics  screens.  All  graphics  are  produced  using  Turbo  Pascal,  a  procedure- 
oriented  computer  language.  In  their  present  form,  the  software  menus  are 
prognmmed  primarily  at  the  DOS  level,  thereby,  making  prog'^am  modifications 
easier.  The  interactive  FORTRAN  programs  automatically  prompt  the  user  tor 
data  input  (on  the  screen)  when  necessary,  and  output  is  sent  both  to  the 
screen  and  to  a  disk  file  that  can  be  readily  printed. 

1.  Binary  Interaction  Databases 

A  “demonstration"  diskette  containing  the  UNIFAC  algorithms  has  been 
developed  to  illustrate  the  potential  of  the  UNIFAC  fragment  approach  for 
calculating  Henry'o  Constants  and  aqueous  solubilities.  The  chemical  fragment 
tables  (Taoles  18  and  19)  accompanying  this  report  should  be  utilized  when 
running  the  UNIFAC  programs.  The  most  recently  published  VLE  and  LlE  data¬ 
bases  are  integrated  into  the  software,  and  either  database  can  be  Inter¬ 
actively  selected  by  the  user  from  the  screen.  UNIFAC  binary  Interaction 
parameters  determined  from  the  expenmenal  Henry's  law  constant  data  have 
also  been  added  to  the  original  VLE  parameters  to  give  a  third  "environmental" 
(infinite  dilution)  database.  The  user,  therefore,  has  three  database  choices 
onscreen  when  running  the  UNIFAC  software. 

Note  from  Tablet  18  and  19  that  the  main  group  and  subgroup  identi¬ 
fication  numbers  for  a  g'ven  fragment  differ  between  the  VLE  and  LLE  data¬ 
bases.  (Subgroup  identification  numbers  for  the  environmental  database  are 
identical  to  those  used  in  the  VLE  database./  The  database  selected  should, 
therefore,  be  kept  in  mind  when  entering  group  identification  nuracers.  In 


TABLE  18.  SUBGROUP  ID  NUMBERS  FOR  VLE  DATA  BASE 


DATABASE  OF  UMIFAC  VLE-PARAHETERS  * 

Group  designationa  of  GaeMlng,  Raaauaaen,  and  Fredenniund,  (19d2>  • 

•  ••  VAPOR-LIQUID  EQULIbRlU!f  •*«  « 


3 


THE  HAIM  GROUPS  ARE; 


1)  CH2 
6)  CH30H 
11)  CCOO  . 
16)  <C)3K 
21 )  CCL  . . 
26)  CH02  . 
31)  DQH  .. 
36)  ACRY  . 


2)  C=C  . 

.  3)  ACH  . . . . . 

.  4)  ACCH2  .. 

. .  5)  OB . 

7)  H20  . 

.  8)  ACOH  .. . . 

.  9)  CP2C0  .. 

.  10)  CKO  .... 

12)  HCOO  .... 

13)  CK20  .... 

14)  CHHE  ... 

IS)  CNN  .... 

17)  ACNH2  ... 

18)  PYRIDINE 

19)  CCN  .... 

20)  COOK  ... 

22)  CCL2  ..  . 

23)  CCL3  .... 

24)  CCL4  ... 

25)  ACCL  ... 

27)  ACHQ2  ... 

28)  CS2 . 

29)  CH3SH  .. 

30)  FURrjRAL 

32)  I . 

33)  BR  . 

34)  C=-C  . . . 

35)  DIISC  '. . . 

37)  CLCC  .... 

38)  ACF  . 

39)  DHF  .... 

40)  CF2  .... 

THE  SUB  GROUPS  ARE; 


1) 

CH3  .... 

...  2) 

CH2 . 

.  3)  CH  . 

.  4) 

c . 

. .  5) 

CH2»CH  ». 

6) 

CH*CH  . . 

. . .  7) 

CH2*C  ... 

.  8)  C4»C  ... 

.  9) 

C»-C  . . . . 

.  18) 

ach . 

11) 

AC . 

.  12) 

ACCH3  ... 

13)  ACCH2  .. 

14) 

ACCH  ... 

.  IS) 

OH . 

16) 

CH30H  .. 

.  17) 

H20 . 

18)  ACOH  ... 

19) 

CH3C0  .. 

.  28) 

CH2C0  ... 

21) 

CHO  .... 

.  22) 

CH3C00  .. 

23)  CK2CCI0  . 

24) 

HCOO  ... 

.  25) 

CH30  .... 

26) 

CH20  ... 

.  27) 

CH-0 _ 

28)  FCH20  .. 

29) 

CH3HH2  . 

.  38) 

CH2NH2  .. 

31) 

CHNH2  .. 

.  32) 

CH3NH  ... 

33)  CH2NH  .. 

34) 

CHNH  ... 

.  35) 

CH3N  .... 

36) 

CK2N  ... 

37) 

ACNH2  .... 

38)  CSHSH  .. 

39) 

CSH4H  .. 

.  40) 

CSH3H  ... 

41) 

CH3CK  .. 

.  42) 

CH2CM  ... 

43)  CCOH  ... 

44) 

HCOOK  .. 

.  45) 

CH2CL  ... 

46) 

CSCl  ... 

.  47) 

CCL . 

48)  CH2CL2  . 

45) 

CHCL2  .. 

.  SC) 

CCL2 _ 

51) 

CKCL3  .. 

.  52) 

CCL3  ....; 

53)  CCL4  ... 

54) 

ACCL  ... 

.  ss; 

CH3N02  .. 

56) 

CH2H02  . 

.  57) 

CHH02  ... 

58)  ACH02  .. 

59) 

CS2  .... 

.  68) 

CK3SH  ... 

SI) 

CH2SH  .. 

.  62) 

FURFURAL 

63)  <CH20H)2 

64) 

I . 

.  65) 

BR . 

66) 

CH—C  .. 

.  67) 

c=-c _ 

68)  DKSO  ... 

69) 

ACRY  ... 

.  78) 

CLiC»-C) 

71) 

ACF  .... 

.  72) 

DKF-1  ... 

73)  DNF-2  .. 

74) 

CF3  .... 

.  75) 

CF2 . 

76)  CF 


[■ABLE  19.  SUBGROUP  ID  NUMBERS  FOR  LLE  DATA  BASE 


•  DATABASE  OF  UHIFAC  LLE-PARAHETERS 

«  Group  deelgnatlons  of  Hognuesen,  Raecaussen,  and  Fredenslund,  (1981) 

•  •••  LIOUID-LIOUID  EOULIBRIUH  ••• 


THE  HAIM  GROUPS  ARE; 


] ) 

rM7  . 

.  7)  r.sr. . 

.  3) 

ach . 

.  4 )  ACCH2  . , 

. . .  5) 

OH . . 

6) 

PI  ... . 

.  . . -  7)  P9 . 

.  8) 

H20 . 

.  9)  ACOH  ... 

. .  10) 

CH2C0  .... 

11) 

CHQ  ... 

..  12)  FURFURAL 

13> 

CCOH  .... 

14)  COOC  ... 

.  15) 

CH20 . 

16) 

CCL  ... 

..  17)  CCL2  ..... 

18) 

CCL3  .... 

19>  CCL4  ... 

.  28) 

ACCL . 

21) 

CCK  ... 

..  22)  ACHH2  ... 

23) 

CH02  .... 

24)  ACN02  .. 

,  25) 

DOH . 

26) 

OEQH  .. 

..  27)  PYRIDME 

28) 

TCE . 

29)  HFA  .... 

.  38) 

D.NFA . 

31) 

THS  ... 

..  32)  DHSO  .... 

THE  SUB  GROUPS  ARE 

1) 

CH3 . 

.  2) 

CH7 . 

.  3) 

CH . 

,  4)  C . . 

.  5) 

CH2«CH  4 

8) 

CH*CH  .... 

.  7) 

CH*C  .... 

.  6) 

CH2*C  . .  - 

.  9)  ach  . 

10) 

AC  ..... 

11) 

ACCH3  .... 

12) 

ACCM2  ... 

13) 

ACCK  ....• 

14)  OH  . 

IS) 

PI . 

16) 

P2 . 

17) 

H20 . 

18) 

ACOH  .... 

19)  CH3CO  ... 

20) 

C82C0  .. 

21) 

CHO . 

22) 

FURFURAL 

23) 

COOH  .... 

24)  HCOOH  ... 

25) 

CK3C00  . 

26) 

CH2C00  ... 

27) 

CH30  .... 

28) 

CH20  .... 

29)  CH-0  _ 

38) 

FCH20  .. 

31) 

CH2CL  .... 

32) 

CHCL  .... 

33) 

CCL . 

34)  CH2CL2  .. 

35) 

CHCL2  .. 

36) 

CCL2 . 

37! 

CHCL3  ... 

38) 

CCL3 _ 

39)  CCL4  .... 

40) 

ACCL  ... 

41) 

CH3CH  .... 

42) 

CH2CH  ... 

43) 

ACKH2  ... 

44)  CH3H02  .. 

45) 

CH2H02  . 

46) 

CHH02  .... 

47) 

ACH02  ... 

43) 

(CH20H)2 

49)  (H0aR)20 

58) 

C5H5N  .. 

SI) 

C5H4K  .... 

52) 

C5H3)!  ... 

53) 

CCL2«CHCL 

54)  HCONHCH3 

55) 

HC0H<CH3)2 

58) 

<CH2)4S02 

57' 

(CH2>2SO 

addition,  enter  only  subgroup  numbers  (not  mam  group  numoers)  and  construct 
the  chemical (s)  of  interest  from  the  appropriate  subgroup  table. 

The  benzene/water  system  is  an  excellent  test  case  for  two  reasons:- 
(!)  that  system's  well-known  experimental  results,  and  (2)  the  simplicity  oT 
subdividing  these  chemicals  into  subgroups,  for  example,  benzene  is  composed 
merely  of  six  aromatic  carbon-hydrogen  subgroups  ("ACH"),  and  water  is  a 
single  subgroup  ("HgO").  Se-/eral  such  test  cases  have  been  run  to  debug  and 
validate  the  UNIFAC  software. 

2.  Documentation  for  the  Simplex  Fitting  Program 

As  described  earlier,  two  FORTRAN  programs  were  written  to  implement 
the  three  fragment  databases:  a  simplex  fitting  algorithm  and  an  interactive, 
menu-dnven  program  to  calculate  Henry's  law  constants  and  aqueous  solubil¬ 
ities.  The  latter  program  can  also  determine  how  a  volatile  chemical  in  a 
multicomponent  liquid  solution  will  partition  between  the  vapor  and  tne  liquid 
phases  at  equilibrium.  Mixture  effects  on  equilibrium  behavior  can,  thus,  be 
quantified  and  predicted. 

The  above  software  items  were  designed  for  ease  of  use.  The  simplex 
routine  was  written  to  fit  UNIFAC  binary  interaction  parameters  to  the  ex- 
per’mental  Henry's  law  constant  data  collected  over  the  course  of  this  study. 
For  the  sake  of  simplicity,  this  program  is  in  "batch"  form  (i.e.,  all  ’nput 
operations  involve  reading  from  ASCII  data  files)  and  is  not  interactive.  The 
simplex  softwarej's  a  modification  of  a  routine  published  by  Fredenslund  et 
a1.  (Reference  1);  the  interested  reader  should  refer  to  their  discussion  of 
the  theory  behind  the  simplex  algorithm.  In  the  procedure  for  running  the 
simplex  program,  the  user  supplies  only  two  things:  (1)  the  subgroup  break¬ 
down  (i.e.,  the  number  of  each  type  of  subgroup)  for  each  chemical  of  in¬ 
terest,  and  (2)  the  experimental  activity  coefficient  Ja.a  for  the  pertinent 
binary  system.  Up  to  four  binary  interaction  parameters  can  be  fitted  simul¬ 
taneously  with  tne  current  software  configuration.  Typically,  these  would  be 
the  interaction  parameters  Aj  ,  and  A,  j  for  groups  "i"  and  "j,"  plus  a  second 
pair,  A|(j  and  Ai_i(,  for  two  additional  subgroups.  Note  that  A,  j  and  Ajj 
are  not,  in  general,  numerically  equivalent.  ' 

As  an  example  of  the  fitting  procedure,  consider  the  determination  of 
the  binary  parameters  for  two  subgroups:  ACH  and  H20.  The  necessary  binary 
experimental  data  could  come  from  any  of  a  number  of  systems  containing  the 
subgroups  of  interest,  such  as  benzene/water,  toluene/water,  or  xylene.-water. 
There  are  a  total  of  three  main  groups  in  these  binary  systems:  (1) 

ACH... (ACH  subgroup!,  (2)  AC-CKZ..,  (AC-CH3  subgroup),  and  (3)  H20...  (H20 
subgroup).  Because  a  main  group’s  interaction  with  itself  is  defined  to  be 
zero,  a  total  of  six  interaction  parameters  must  either  be  fitted  or  supplied 
by  the  user.  The  toluene/water  system,  tor  instance,  has  every  one  of  the  six 
(6)  possible  binary  interaction  pairs- 

•  ACH  and  AC-CH3, 

•  AC-fH3  and  ACH, 
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ACH  and  H20 


•  H20  and  ACH, 

•  AC-CH3  and  H20, 

•  H20  ana  AC-CH3. 

As  further  illustration,  the  simplex  algorithm  could  simultaneously 
fit  interaction  Parameters  3-6  above,  but  values  for  Parameter  1  and  2  are 
needed  to  complete  the  structural  description  of  both  toluene  and  xylene. 
Parameters  1  and  2  would  either  be  estimated,  supplied  from  a  previous  simplex 
fitting  run,  or  taken  from  the  published  VLE  database  of  Gmehling,  Rasmussen, 
and  Fredenslund  (Reference  9).,  Clearly,  when  multiple  compounds  and  subgroups 
are  considered,  the  fitting  procedure  becomes  very  complex.  In  the  present 
work,  a  sequential  fit  (i.e,  using  previously  fitted  interaction  pairs  to  a’d 
in  fitting  subsequent  pairs)  was  selected  to  minimize  computer  time.  Perhaps 
a  more  desirable  future  approach  would  be  to  fit  all  interaction  parameters 
simultaneously  with  a  pattern  search  algorithm  (e.g.,  Hooke-Jeeves)  to  mini¬ 
mize  the  difference  between  the  experimental  and  calculated  activity 
coefficients. 

3.  Documentation  tor  the  Interactive  UNIFAC  Program 

The  interactive  UNIFAC  program  to  calculate  Henry's  law  constant  and 
aqueous  solubility  is  very  simple  to  operate.  It  can  be  run  on  any  IBM-com¬ 
patible  personal  computer  that  supports  the  MS-DOS  or  PC-DOS  operating  systems 
and  may  be  installed  on  a  hard  d'sk  if  desired.  To  run  the  program  from  a 
floppy  disk,  the  following  steps  should  be  performed: 

•  Insert  the  program  disk  into  the  designated  drive  and  change  the 
screen  drive  prompt  to  the  appropriate  letter  (e.g.,  "A>"). 

•  Type  the  command  “UNIFAC  at  the  DOS  drive  o'ompt  and  press  [Enter] 
to  invoke  the  program. 

•  A  graohics  start-up  screen  will  apoear  while  the  program 
initializes,  after  which  a  selection  menu  will  appear. 

•  Select  "A"  for  calculation  of  Henry's  Constant,  "8"  for  the  simplex 
routine  (not  installed  in  the  current  version),  and  "C’  for  the 
liquid-liquid  flash  algorithm  (e.g..  an  iterative  Newton-Raphsoii 
procedure)  to  calculate  aqueous  solubility. 

•  Follow  all  instructions  as  they  are  displayed  onscreen  and  supply 
information  from  the  keyboard  as  requested.  Program  output  is 
written  both  to  the  console  (screen)  and  to  a  disk  ASCII  file  that 
can  be  routed  to  virtually  any  printer  or  imported  into  Lotus  1-2-3 
and  most  other  commercial  programs.  Printing  the  program  outout  can 
be  done  with  the  following  DOS  copy  command:  "copy  UNIFAC. OUT 
IptI:."  This  assumes  connect-on  o^  the  printer  to  parallel  port  1, 
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but  other  printer  interfaces  (e.g.,  Ipt2:,  coml:,  and  corri2:)  can 
also  be  used. 

•  After  calculations  are  terminated  and  the  results  displayed,  the 
program  queries  the  user  whether  to  return  to  the  Main  Menu  or 
continue.  This  feature  allows  the  user  to  re-run  a  branch  of  the 
program  multiple  times  without  ha/ing  to  return  to  the  Ham  Menu 
after  each  pass. 

•  At  the  Ham  Menu,  the  user  has  the  option  of  selecting  another  item 
or  exiting  the  program  by  pressing  the  [Escape]  key.  Leaving  the 
program  returns  tne  user  to  the  initial  DOS  orive  prompt. 

To  operate  the  program  from  a  hard  disk  instead  of  a  floppy  disk, 
simply  create  a  hard  disk  directory  (at  the  "C>"  prompt)  with  the  DOS  command 
"md  c:\UNIFAC.''  Then  place  tne  original  program  diskette  into  drive  A  and 
enter  the  following  DOS  coiimand  at  the  "O"  prompt:  "copy  a:*.*  c:\L'NlFAC." 

To  run  the  program,  simply  enter  the  DOS  command  "cd  \LINIFAC"  at  the  "C>" 
prompt  and  type  "UNIFAC."  The  program  should  tnen  run  normally  according  to 
the  above  instructions  for  floppy  disk  operation.  When  duplicating  the 
original  program  diskette,  always  copy  the  contents  of  the  entire  diskette 
since  the  UNIFAC  program  uses  all  the  files  and  expects  to  find  them  m  the 
startup  directory.  Because  the  original  UNIFAC  diskette  contains  the  MS-DOS 
2.11  operating  system,  it  is  recommended  that  the  "diskcopy"  commands  be  used 
to  exactly  duplicate  the  UNIFAC  disk.  The  original  program  diskette  will  then 
"boot'  the  system.  That  is,  if  the  disk  is  placed  in  drive  A:  and  the  com¬ 
puter  turned  on,  the  UNIFAC  program  will  come  up  without  any  further  commands. 
For  users  with  limited  computer  experience,  this  may  be  the  preferred  call-up 
procedure. 

The  only  data  necessary  to  run  the  interactive  program  are  the  sub¬ 
group  ID  numbers  (see  Table  18  or  19)  that  describe  the  chemicals  of  interest, 
the  absolute  temperature  of  the  system,  and  pure  component  vapor  pressure  data 
if  Henry's  law  constants  are  to  be  calculated.  Antoine  coefficients  can  also 
be  supplied  instead  of  component  vapor  pressures  if  they  are  available. 


SECTION  VII 


CONCLUSIONS 


Experimental 1y  determined  Henry's  law  constants  from  this  study  agreed 
well  with  other  reported  values.  Howe- er,  for  many  of  the  compounds  there 
were  very  few,  if  any,  experimentally  determined  values  with  which  to  compare, 
and  rarely  were  the  results  presented  with  confirming  statistics.  In  many 
cases,  agreement  between  experimentally  determined  Henry's  law  constants  and 
values  estimated  using  the  ratio  of  vapor  pressure  and  aqueous  sol ubi 1 ity  was 
quite  good.  However,  depending  on  the  application  of  the  data  (i.e.,  required 
accuracy)  car^  should  be  exercised  when  using  estimated  values,  as  differences 
of  over  400  percent  were  noted.  Equilibrium  Partitioning  in  Closed  Systems 
(EPICS)  is  an  attractive  methodology  for  directly  determining  these  part'tion- 
ing  values.  This  simple  technique  employs  standard  gas  chromatography  head- 
space  analysis,  and  absolute  gas  phase  concentrations  are  not  needed  for  the 
analysis.  The  EPICS  headspace  method  is  applicable  over  a  wide  volatility 
range  and  exhibits  excellent  precision  and  reproducibility  from  one  laboratory 
to  another.  Batch  air  stripping  is  comparable  in  accuracy  to  EPICS  for  dilute 
aqueous  solutions  of  less  volatile  compounds.  However,  long  bubble  retention 
times  may  be  required  to  achieve  gas-liquid  equilibrium  for  volatile  constitu¬ 
ents,  and  axial  concentration  gradients  may  arise  when  tall  columns  or  low  gas 
flow  rates  are  employed. 

Mixtures  of  different  contaminants  may  cause  both  positive  and  negative 
deviations  from  ideal  solution  behavior  due  to  intermolecular  effects  such  as 
association  and  solvation.  These  systems  of  mixtures  more  closely  resemble 
actual  environmental  contamination  problems  than  do  neat  single  component 
solutions.  The  EPICS  technique  allows  ready  analysis  of  the  apparent  Henry's 
law  constant  of  each  component  in  these  mixtures.  This  information  would 
permit  the  more  reliable  design  of  treatment  facilities  for  contaminated 
groundwater  or  other  waste  matrices;  however,  meaningful  interpretation  of 
these  results  awaits  the  development  of  a  predictive  multicomponent  model. 

The  shake-flask  method  proved  to  be  the  most  reliable  technique  for  deter¬ 
mining  aqueous  solubilities  of  the  three  methods  employed.  Local  maxima  and 
minima  of  solubility  with  respect  to  temperature  were  observed,  but  often  these 
extremes  were  not  statistically  significant.  The  correlation  of  the  solubility 
data  to  the  Henry's  law  constant  results  is  not  possible  at  this  time. 

The  UNI  FAC  method  has  proven  to  be  the  most  effective  way  of  utilizing  the 
current  VLE  data  in  a  general  thermodynamic  correlation  of  aqueous  solubilities 
and  Henry's  law  constants.  A  computer  algorithm  to  fit  the  current  data  to  a 
new  environmental  UNIFAC  binary  interaction  database  has  been  developed,  and  a 
portion  of  the  experimental  data  collected  incorporated  into  this  new  database. 
The  new  database  creates  a  considerable  improvement  in  the  predictions  gener¬ 
ated  by  UNIFAC  in  the  dilute  concentration  regime  (Figures  13  through  16). 
Additional  effort  will  be  required  to  completely  fit  the  current  experimental 
data  with  the  environmental  UNIFAC  database. 
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SECTION  VIII 


RECOMMENDATIONS 


Based  on  the  encouraging  results  of  this  study,  there  are  opportunities 
for  future  work  in  several  areas:  (1)  extension  and  refinement  of  the  binary 
interaction  fitting  procedure  begun  in  this  study,  (2)  continued  development 
of  interactive  software  to  implement  the  UNIFAC  concept,  and  (3)  additional 
experimental  measurements.  Continued  work  in  all  three  areas  is  essential  to 
achieving  the  goal  of  a  comprehensive  UNIFAC  "environmental"  database  that  is 
easy  to  use. 

A.  CORRELATION  DEVELOPMENT 

Further  work  in  fitting  UNIFAC  binary  interaction  parameters  would  be 
greatly  aided  by  the  following  actions: 

•  Extension  of  the  simplex  program  to  handle  the  simultaneous 
fitting  of  more  than  the  current  two  pairs  of  interaction 
parameters 

•  Modification  of  the  simplex  program  to  eliminate  forced 
optimization  in  pairs,  especially  if  the  water-subgroup 
interaction  parameter  (e.g.,  Ah20-***)  is  assumed  to  be 
invariant 

•  Rewrite  the  simplex  program's  input-output  statements  to 
eliminate  the  necessity  of  entering  water  mole  fractions  and 
experimental  activity  coefficients  since  these  are  the  same  in 
all  cases 

•  Change  the  output  file  format  for  the  interactive  Henry's  law 
program  to  match  the  necessary  input  file  format  for  the  simplex 
program,  and  perhaps  make  this  common  format  compatible  with  a 
standard  plot  routine. 

Computer  software  development  is  an  essential  part  of  implementing  the 
UNIFAC  algorithm  in  a  manner  that  meets  Air  Force  needs.  The  software  package 
developed  under  the  current  contract  can  be  improved  in  several  specific  ways. 
One  area  of  future  attention  would  be  better  screen  cursor  positioning  and 
data  input.  For  instance,  the  program  now  prompts  the  user  sequentially  for 
subgroup  identification  numbers  and  does  not  allow  an  entered  value  to  be 
changed.  A  better  approach  would  be  to  have  full-screen  input  of  data  by 
columns. 

Applied  to  Subgroup  ID  Inputs,  this  would  Involve  displaying  the  available 
subgroups  on  the  screen  in  several  parallel  columns  with  a  highlighted  cursor- 
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arrow  that  can  be  positioned  by  the  user  next  to  any  desired  group.  The  num¬ 
ber  of  each  type  of  subgroup  in  a  given  molecule  could  then  be  typed  next  to 
the  appropriate  subgroup  ID  on  the  screen.  A  menu  selection  for  "Resume  Pro¬ 
gram  Execution"  would  be  selected  by  the  user  only  after  the  subgroup  break¬ 
down  of  the  molecule  is  completely  defined..  While  on  the  screen  (befo'-e  final 
selection),  an  entry  could  be  changed  at  any  time  by  the  user.  These  and 
other  useJful  improvements  to  the  interactive  software  could  not  be  implemented 
within  the  scope  of  the  present  contract. 

S.  FUTURE  EXPERIMENTAL  WORK 

Because  pure  component  vapor  pressure  data  are  necessary  to  calculate 
Henry's  law  constants  from  UNIFAC  activity  coefficients,  additional  experi¬ 
mental  work  should  include  the  measurement  of  saturation  pressures  as  neces¬ 
sary.  For  those  compounds  whose  published  Antoine  constants  do  not  apply  at 
near-ambient  environmental"  temperatures  (or  are  not  available),  at  least  one 
saturation  pressure  measurement  between  10  'C  and  30  *C  is  recommended. 

Experimental  determination  of  Henry's  law  constants  and  aqueous  solubil¬ 
ities  should  oe  designed  to  complement  the  thermodynamic  modeling  objectives. 
For  instance,  missing  UNIFAC  interaction  parameters  between  subgroups  can  be 
used  to  logically  select  experiments  that  fill  such  date  gaps.  For  this  rea¬ 
son,  future  experiments  should  focus  on  oxygenated  compounds  and  chemicals 
with  multiple  ponds  (e.g.,  olefins).  Such  chemicals  would  be  environmentally 
interesting  in  their  own  right  and  would  increase  the  scope  of  the  environ¬ 
mental  UNIFAC  database  being  developed. 
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